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We review several cosmological backgrounds of gravitational waves accessible to direct-detection 
experiments, with a special emphasis on those backgrounds due to first-order phase transitions 
and networks of cosmic (super-)strings. For these two particular sources, we revisit in detail the 
computation of the gravitational wave background and improve the results of previous works in the 
literature. We apply our results to identify the scientific potential of the NGO/eLISA mission of 
ESA regarding the detectability of cosmological backgrounds. 
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I. INTRODUCTION 



Gravitational waves represent one of the most promising messengers to probe the Universe in some of its most 
fundamental manifestations, in particular the violent phenomena associated with the last stages of the evolution 
of stars, or with the collisions of galaxies and the early history of the Universe. Various frequency windows of 
gravitational radiation are now being targeted by different experimental means, in particular [10 _18 Hz, 10~ 16 Hz] 
through the polarization of the Cosmic Microwave Background, [10 _9 Hz, 10 _7 Hz] through millisecond pulsar timing, 
[10 -4 Hz, 10 _1 Hz] with space interferometers and [10Hz, 10 3 Hz] with ground interferometers. The LISA program has 
been pursuing the goal of detecting gravitational waves in the sub-Hz region for more than ten years as a joint venture 
between ESA and NASA. A new development in this program was the decision by Europe early 2011 to move forward 
on its own. The new Europe-only mission, named New Gravitational wave Observatory (NGO) by ESA and eLISA 
(evolved LISA) by the community 1 , departs significantly from the previous ES A-NAS A mission. It is thus important 
to assess its scientific program. This important endeavour has been undertaken by various groups in the last months, 
covering the full spectrum of the scientific program of the new mission |l| . We report here on the assessment of the 
scientific potential of the new mission regarding the measurement of cosmological backgrounds. Although not the 
main motivation for the mission, the detection of such backgrounds would be a major discovery. Indeed, gravitational 
waves (GW) penetrate all of cosmic history, which allows GW detectors to explore scales, epochs, and new physical 
effects not accessible in any other way [2|-l5|- Two important mechanisms for generating stochastic backgrounds are 
phase transitions in the early Universe and cosmic strings. 

Gravitational waves produced after the Big Bang form a fossil radiation: expansion prevents them from coming in 
thermal equilibrium with the other components because of the weakness of the gravitational interaction. Important 
information on the first instants of the Universe is thus imprinted in these relics and can be decoded. The mechanical 
effect of expansion is simply to redshift the corresponding frequency (and amplitude). Assuming that the wavelength 
is set by the apparent horizon size c/H* = c(a/a)* at the time of production (when the temperature of the Universe 
is T*), the redshifted frequency today is 



Thus, eLISA frequency band of about 0.1 mHz to 100 mHz today corresponds to the horizon at and beyond the 
Terascale frontier of fundamental physics. This allows eLISA to probe bulk motions at times about [3 • 10 -18 , 3 • 10~ 10 ] 
seconds after the Big Bang, a period not directly accessible with any other technique. Taking a typical broad spectrum 
into account, eLISA has the sensitivity to detect cosmological backgrounds caused by new physics active in the range 
of energy from 0.1 TeV to 1000 TeV, if more than a modest fraction of about 10~ 5 of the total energy density is 
converted to gravitational radiation at the time of production. A standard example of new physics is a first-order 
phase transition resulting in bubble nucleation and growth, and subsequent bubble collisions and turbulence. 

Phase transitions also often lead to the formation of one-dimensional topological defects known as cosmic strings. 
Among possible topological defects, cosmic strings are unique from a cosmological point of view because, whereas their 
relative energy density should grow with the expansion, they interact and form loops which decay into gravitational 
waves. Thus cosmic strings tend to form networks with a typical scaling behaviour, losing energy mainly through 
gravitational radiation with a very broad and uniquely identifiable spectrum. Besides topological defects in field 
theories, cosmic strings may also occur as fundamental objects of string theory, the theory that is aiming to provide 
a unified framework for all particles and forces of nature. Indeed, although fundamental strings were devised as sub- 
microscopic objects, it has been progressively realized [y] that some of these strings could be stretched to astronomical 
size by the cosmic expansion. eLISA will be a very sensitive probe for these objects and so offers the possibility of 
detecting direct evidence of fundamental strings. 

We review in what follows the different cosmological backgrounds that eLISA could potentially probe, with a special 
attention to first-order phase transitions and cosmic strings. Because the study of these backgrounds has made a lot 
of progress over the years, we review the literature in a rather detailed way in order to motivate our own choices of 
parameters and dynamical processes at work, choices which may not coincide with previous analyses. The aim of this 
work is to provide reliable and precise detection forecasts for the GW signal coming from these two sources. In order 
to do that, we have improved and completed the models of the GW spectra available in the literature in the following 
aspects. 

In the case of first order phase transitions, the GW signal has been evaluated in previous analyses under the 
assumption that the broken phase bubbles propagate through a Jouguet detonation. This assumption simplifies 



In what follows, we use the latter name. 
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considerably the evaluation of the signal but is very restrictive and not justified in the vast majority of phase transition 
scenarios. Here this approximation is relaxed: as a consequence, we find that the amplitude of the GW signal depends 
not only on the temperature, strength, and duration of the phase transition, but also on a new parameter representing 
the friction exerted on the bubble walls by the interaction with the surrounding plasma. Moreover, we combine for 
the first time the GW spectra from bubble collisions and magnetohydrodynamic turbulence in a fully consistent way, 
and provide detection forecasts for the total GW signal. In previous works on the estimation of the GW signal from 
first order phase transitions, the strength and the duration of the phase transition have often been considered as 
independent parameters. However, a relationship among them is to be generically expected for thermal first-order 
phase transitions: we account for this fact in the present evaluation of the GW signal, and point out that it has 
relevant consequences as far as detection prospects are concerned. 

In the case of cosmic (super-)strings, two main methods have been used in the literature to evaluate the GW 
background. The first one is based on models for the GW spectrum emitted by each individual loop, while in the 
second method the background is obtained as a superposition of GW bursts produced by specific loop configurations 
called cusps and kinks. We show that the two methods lead to similar results up to two differences: an overal 
normalisation of the GW background and the effect of removing the rare bursts that do not contribute to the gaussian 
and stationary background. We study carefully how these two differences affect the results. We also show that the 
GW background today depends weakly on the particular spectrum emitted by each individual loop. For instance, 
contrary to what is often stated in the literature, cusps do not lead to any significant increase of the background 
amplitude. On the other hand, the signal depends strongly on properties of the string network evolution that are still 
uncertain, in particular the characteristic size of the loops when they are produced. We therefore study in detail how 
these uncertainties affect the results. Compared to previous works, we also improve the modelling of the cosmological 
evolution and we study how the results depends on the thermal history of the very early universe. Indeed, since 
cosmic strings emit GW continually during the cosmological history, the GW background today depends also on the 
details of the cosmological evolution. 

Readers interested only in the scientific potential of eLISA may go directly to Section III D 21 for first-order phase 
transitions and to Sections lIII C 3l and lIII D 3l for cosmic strings. For a brief but exhaustive review of other cosmological 
sources of GW, see Section ITVl 

The rest of the paper is organized as follows. We briefly review general aspects of cosmological GW backgrounds 
in sub-section II Al introducing our notations, and discuss their detection with eLISA in sub-section II Bl Section [H] 
is dedicated to first-order phase transitions. After an general overview of the GW sources and spectrum (sub- 
section III Ap , we review the literature about the calculation of the GW background (sub-section III Bp and about 
particle-physics models of first-order phase transitions at the electroweak scale ( sub-section III C[) . In sub-section III Dl 
we then use the latest results available in the literature to evaluate the GW spectra in two models of first-order 
phase transition, and to assess more generally the detection prospects for eLISA. Section IIIII is dedicated to cosmic 
(super-)strings. We start in sub-section 1111 Al with an introduction and overview of the literature, where we review 
the main aspects that must be modeled in order to compute the GW background. In sub-section IIII Bl we detail 
our own method to calculate the GW spectrum and compare it to previous methods used in the literature. The 
cases of small and large initial loop sizes are then discussed separately in sub-sections IIII CI and IIII Dl where we 
check the dependence of the results on the model assumptions and determine the regions of the parameters space 
that are accessible to eLISA, as well as ground-based experiments and pulsar timing observations. In Section [IVI we 
review several other cosmological sources of GW, that operate during inflation (sub-section IIV A| . just after inflation 
(sub-section IIV B"|) and during the subsequent thermal evolution (sub-section IIV Cj) . Finally, Section IVl contains our 
conclusions. 

A. General Aspects and Notations 

In a cosmological context, GW may be represented by a tensor perturbation hij = 1,2,3) of the Friedmann 
Robertson- Walker metric 

ds 2 = -dt 2 + a 2 {t) (8ij + 1 1,, ) d.r'd.r' (2) 

which is transverse and traceless 

dihij =hu = 0, (3) 

where we assume flat spatial sections, a(t) is the scale factor, t denotes the physical time and repeated latin indices 
are summed. The transverse-traceless condition ([3]) leaves only two independent degrees of freedom, which are the 
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only ones that propagate and carry energy out of a source. In Fourier space, their linearized equation of motion given 
by Einstein equations is 2 

hij (k, t) + 3H ha (k, t) + — h i:i (k, t) = 16ttG n|J T) (k, t) , (4) 

(TT) 

where G is Newton constant, H is the Hubble rate, a dot denotes derivative with respect to t, and n^- is the 
transverse-traceless part of the anisotropic stress . The latter is given by a 2 Hij = Tij — pa 2 (5ij + hij ) , where Tij 
denotes the spatial components of the energy-momentum tensor and p is the background pressure. In Eq. ((4]) and in 
the following, we use k = |k| to denote the comoving wave-number of the GW, while we will use k — k/ato denote the 
physical wave-number at a given moment of time. For a source operating at sub-Hubble scales (k 3> H), Eq. ((4]) may 

be approximated as a standard wave-equation 3 hij + k 2 hij = \QttG T^J T ^ / a 2 . Another mechanism of GW production 
is provided by the parametric amplification of quantum fluctuations during inflation. In that case, the source term is 
negligible in Eq. (j4]), but super-Hubble tensor perturbations are generated by the very fast expansion of the universe. 
These tensor perturbations become standard GW once they re-enter the Hubble radius during the post-inflationary 
evolution. 

Once produced in the early universe, GW propagate essentially freely until today, being simply redshifted by the 
expansion of the universe. Their energy density today can be written as [3| 

(hij hij) _ f df dp gw 



Pgw ~ ~~32^G~ ~ J J dh^f ' {b) 

where / = k/(2irao) is the present-day GW frequency (do denoting the scale factor today) and () denotes ensemble 
average. The superposition of GW produced by a large number of unresolved sources in the early universe form a 
stochastic background that is assumed to be statistically isotropic, stationary and nearly Gaussian 2]. Its main prop- 
erties arc then described by its power spectrum. The quantity that is usually considered to characterize cosmological 
backgrounds is the spectrum of energy density per logarithmic frequency interval divided by the critical density p c 
today 

h 2n gw (f)=(^^-) (6) 

\Pc dlogf J Q 

where the subscript refers to the present epoch and h parametrizes the small uncertainty in the value of the Hubble 
constant today, H = 100 /ikm/s/Mpc. 

A source that operates at sub-Hubble scales at some time £* after inflation emits GW with a characteristic physical 
wave- number fc* = k/a* that is larger than the Hubble rate if* at that time: fc* = ifj/e, with e* < 1. The 
characteristic GW frequency today is then given by f c = (k*/2n) (a* /do). For GW produced in the radiation era 
when the plasma temperature is T*, and assuming a standard adiabatic thermal history for the evolution of the 
universe after GW production, the characteristic frequency today can be written as 

/t ^ft)[£)'" (r) 



^lTeVy VlOO/ 

where is the number of relativistic degrees of freedom at temperature T*. The parameter e* < 1 depends on 
the dynamics of the particular GW source under consideration. For a first-order phase transition, one may have for 
instance e* ~ 10~ 3 — 1. In that case, Eq. ([7]) shows that GW produced around the electro- weak scale are potentially 
interesting for detection with eLISA. On the other hand, cosmic string loops produce GW continually (over a wide 
range of values of T*), so that the present-day spectrum covers a very wide range of frequencies. For a given frequency 
of observation, the signal is dominated by the emission at some time i* from loops of size smaller than the Hubble 
radius by many orders of magnitude, so the typical values of e* are much smaller in that case. As another example, 
Eq. ([7]) can also be applied to GW from inflation, but with T* the temperature when a given tensor mode re-enters 
the Hubble radius after inflation, hence with e* = 1. The resulting GW spectrum covers also a very wide frequency 
range, since tensor modes continually re-enter the Hubble radius. 



2 From now on, we set kg = c = 1. 

3 The terms in hij and 16nGphij in Eq. l[4j are negligible compared to k 2 hij for A: 2> H. 
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FIG. 1: Sensitivity curve of eLISA computed using the expected instrumental noise and the confusion noise generated by 
unresolved galactic binaries [lj. 

Before moving on to eLISA, let us briefly review current constraints on cosmological GW backgrounds in different 
frequency ranges, see e.g. [j| and references therein. If GW contributed too much to the total energy density of the 
universe at the epoch of Big Bang Nucleosynthesis (BBN) , when T ~ MeV, they would spoil the succesfull predictions 
for the light element abundances. This leads to a constraint on any GW background that is present at that time, 
which can be written in terms of the present-day GW energy density as h 2 (p gw /Pc)o < 5.6 x 10 -6 (N v — 3), where N v 
is the effective number of neutrino species. Although this bound applies to the total energy density in GW, neglecting 
the unlikely possibility of a GW spectrum with a very narrow peak one can take as a rule of thumb h 2 Q, gw < 10~ s 
for / ;> 1CP 10 Hz (/ ~ 10~ 10 Hz corresponds to the comoving Hubble scale at the time of BBN). Cosmological GW 
at very low frequencies are also constrained from their contribution to the temperature anisotropies of the Cosmic 
Microwave Background (CMB). This leads to the so-called COBE bound for GW that entered the Hubble radius after 
the time of last scattering, corresponding to frequencies in the range / ~ 1CP 18 — 10~ 16 Hz. The strongest bound 
on h 2 Sl gw is obtained at / ~ 10~ 16 Hz: h 2 Vt gw < few x 10~ 14 . A GW background with a too large amplitude can 
also perturb the very accurate timing of msec pulsar. The resulting constraint obtained in 0] is h 2 VL gw < 2 x 10~ 8 
at / = l/(8yr) (up to a weak dependence on the slope of the spectrum). Finally, the observational upper bound 
h 2 Q. gw < 3.6 x 1(T 6 at / w 100 Hz has been obtained from the S5 LIGO run 0. 

B. Looking for Stochastic Backgrounds with eLISA 

NGO/eLISA consists, like the original LISA mission, of a triangular constellation of 3 satellites. But whereas LISA 
had 6 laser links joining the three satellites (two in opposite directions for each side of the triangle) , eLISA only has 
four lasers, which link one satellite (called Mother) to the other two (called the Daughters). Hence, eLISA consists of 
a single Michelson-type experiment. The inter-spacecraft distance is reduced to one million kilometers (down from 5 
millions in the case of LISA) and the nominal mission lifetime to 2 years (extendable to 5 years). The eLISA sensitvity 
curve in terms of h 2 fl gw is given in Figure [TJ 

In order to distinguish backgrounds of gravitational waves from waves emitted by point sources, it is essential 
to make use of the successive positions of eLISA around the Sun, and thus to wait a sufficient amount of time (of 
the order of a few months). Moreover, it will be more difficult for eLISA, with respect to LISA, to disentangle 
an isotropic cosmological (or astrophysical) background from an instrumental one, all the more because the eLISA 
Mother-Daughter configuration, providing only two measurement arms, does not allow to use Sagnac calibration 
Luckily as we will see, in the case of phase transitions as well as cosmic strings, the spectral dependence of the signal is 
well predicted and may allow to distinguish cosmological backgrounds as long as they lie above the eLISA sensitivity 
curve. 
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II. FIRST-ORDER PHASE TRANSITIONS 
A. General Overview: Characteristic Frequency and Amplitude of the GW Signal 

In the course of its adiabatic expansion, the universe might have undergone several phase transitions (PTs) driven 
by the temperature decrease. The nature of the primordial PTs depends on the particle theory model, but if they are 
first-order they proceed through the nucleation of broken phase bubbles, which is a very violent and inhomogeneous 
process capable of sourcing GW. In this section we concentrate specifically on GW generation by processes which are 
related to bubble nucleation during a first-order PT. We study the characteristics of the GW power spectrum and its 
dependence on the PT parameters, and provide a complete analysis of the GW signal improving the results present 
in the literature, in particular concerning bubble propagation and the energy balance of the PT. 

The electro- weak phase transition (EWPT) in the standard model is a crossover, and it is not expected to lead to any 
appreciable cosmological signal; however, deviations from the standard model in the Higgs sector or the introduction 
of supersymmetry can lead to a first-order EWPT. Similarly, the QCDPT is also predicted to be a crossover by lattice 
simulations but it can become first-order if the neutrino chemical potential is sufficiently large [10| . GW detection 
would help to probe the nature of these PTs, and provide interesting information on the underlying particle theory. 
In particular, a sufficiently strong first-order EWPT is potentially very interesting, because it could have produced a 
relic GW signal right in the frequency band of eLISA (c.f. section IT"A"|) . 

Towards the end of the PT, the true vacuum bubbles collide and convert the entire universe to the broken phase. 
The collisions break the spherical symmetry of the bubble walls, generating a non-zero anisotropic stress which acts 
as a source of GW (see Eq. Moreover, bubble collision causes an injection of energy in the primordial plasma, 

which has a very high Reynolds number (of the order of 10 13 at 100 GeV and at the typical scale of the bubbles 
fill]): the energy injection leads to the formation of magnetohydrodynamic (MHD) turbulence, which sources GW 
through the anisotropic stresses of the chaotic fluid motions. MHD turbulence also leads to the amplification of small 
magnetic fields generated by charge separation at the bubble wall [l2l . fl3j j , which have non-zero anisotropic stress and 
also source GW. 

There are in summary two processes which can lead to the production of GW towards the end of a first-order PT: 
bubble collision and MHD turbulence (where the anisotropic stress is due both to the velocity field and to the magnetic 
field). They are related to the collision of bubbles, which involves two quantities: the duration of the PT, commonly 
denoted by the parameter and the typical size of the bubbles at the moment of collision, i?* ~ Vb/3" 1 , where Vb 
is the bubble wall velocity. The characteristic frequency of the GW generated by the two processes can correspond 
either to the duration of the PT or to the bubble size: fc* ~ f3 or i?" 1 , depending on the details of the time evolution 
of the source (c.f. section Hi B|) . These two parameters can therefore both determine the value of e* in Eq. (J7)). If 
the growth of the bubble proceeds at a highly relativistic speed, the two time/length-scales are equal. Assuming for 
the moment fc* ~ /?, one obtains the following order of magnitude estimate of the characteristic frequency today (c.f. 
Eq. 0): 

The parameter /3/H* is the ratio of the Hubble time to the duration of the PT. Since the entire universe must be 
converted to the broken phase, the PT must complete faster than a Hubble time, so in general we expect /3/H* > 1. 
From Eq. (j8]) it appears that the characteristic frequency of GW emitted at the EWPT at 100 GeV falls in the 
frequency range of eLISA for values 1 < /3/H* < 10 4 . As another example, we see from the above formula that GW 
production at the QCDPT at T» ~ 100 MeV can fall into the frequency range of detection with pulsar timing array, 
/ > 10~ 8 Hz (see e.g. [3] and references therein). The precise value of /3/-H* has to be determined in the context of 
a given model for the first-order PT (c.f. section Hi Dp . 

A simple estimate of the GW amplitude, which shows how the result scales with the duration and the energy 
density of the GW source, can be given through the following heuristic argument. We rewrite the GW equation of 
motion (|4]) for > 1 simply as f3 2 h ~ WirGT, where h denotes the amplitude of the tensor perturbation, T the 

energy-momentum tensor of the source, and we inserted 1//3 as the characteristic time on which the perturbation is 
evolving (we have dropped indices for simplicity). This suggests that h ~ 16irGT/f3, and the GW energy density 
at the time of production can then be estimated as p gw ~ h 2 / (32irG) ~ 8TrGT 2 /f3 2 . Dividing by the total energy 
density p tot — 3H 2 / (8ttG) at the time of GW production, we can write p gw j ' p tot ~ (-ff*//3) 2 (p s /ptot)t where p s ~ T 
denotes the part of the energy density available in the source for the GW generation. Accounting for the fact that 
the PT takes place in the radiation-dominated universe and that the GW energy density is diluted like radiation, one 
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then obtains, for the peak amplitude of the GW spectrum (J6j) today 

where £Ir = pr/ p c denotes the radiation abundance today and {p s /ptot)* the fraction of energy density that con- 
tributes to GW generation at the time of production. The above equation shows that the GW energy density scales 
like the square of the ratio of the GW source duration and the Hubble time, and the square of the energy density 
in the source. As a rule of thumb, given h 2 flji ~ 4.15 x 10~ 5 , a GW signal above the lowest sensitivity of eLISA 
(h 2 Q gw > 10~ 10 ) can be generated if (if*/ '/3)(p s / ptot)* ^ 3 x 1CP 3 . Therefore, detectable signals arise from very 
energetic processes, which involve a sizable fraction of the total energy density in the universe, and at the same time 
slow processes, which minimise the value of fi/H* (in compatibility with the value of the characteristic frequency, 
Eq. which should not exit the detection frequency band - c.f. the discussion in section Hi D I and the spectra in 
Fig.©. 

The slope of the GW spectrum at wave-numbers fc smaller than the Hubble radius H* at the time of production can 
also be determined on general grounds, valid for any transient stochastic source after inflation. This is a consequence 
of the fact that the causal process generating the GW signal cannot operate on time/length-scales larger than l/fi„. 
Therefore, the anisotropic stresses TL-J (k, f) sourcing the metric perturbations in Eq. (|4| are not correlated for 
k < if», and the anisotropic stress power spectrum is expected to be flat (white noise) up to the wavenumber fc* 
(beyond fc*, the power spectrum can instead decay with a slope that depends on the details of the source, see e.g. [l5| 
and references therein). If the anisotropic stress power spectrum is flat, so is expected to be the GW power spectrum 
\hij(k, t)\ 2 , because the causal horizon sets also the maximal time-scale of the correlation, and Eq. (|4]) cannot lead 
to any extra correlation for k < H*. This in turn implies that the spectrum of GW energy density per logarithmic 
frequency interval must grow as fc 3 , because dp gw /d\ogk cx k 3 \hij(h,t)\ 2 . Thus the infra-red tail of the present-day 
GW spectrum behaves as h 2 £l gw cx f 3 for scales that were super- Hubble at the time of production. 

The rest of this section|TT]on GW generation by first-order PTs is organised as follows: we start with an overview of 
the literature, in order to present the main achievements and resume the state of the art of the subject (sections III Bl 
and III Cl) . We then explain how the results of the literature can be improved, in particular concerning two aspects: 
the shape of the GW spectrum arising from MHD turbulence, and the estimation of the PT parameters, especially 
about the physics of bubble propagation (section |H Dp . We point out that restrictive formulas have been used for the 
bubble wall velocity and the kinetic energy of the bubble wall, and use new results of the literature to correct this 
problem (section III D II) . In the last section, HID 21 we evaluate the GW spectra in two examples of first-order PTs, 
and then discuss the prospects for detection with eLISA in a way as mo del- independent as possible. 

B. Overview of the Literature: GW Production by Bubble Collision and MHD Turbulence 

1. Bubble Collision 

The fact that a first-order PT occurring explosively, throu gh the nucleation of fast broken phase bubbles, can be 
a source of GW, has been first pointed out in Witten 1984 This work estimates the characteristic frequency 

and intensity of GW emitted during a first-order QCDPT, finding that the frequency falls in the sensitivity range of 
pulsar timing array observations: if i?* denotes the radius of the typical broken phase bubbles, [l6| demonstrates for 
the first time that the wavelength of the emitted GW corresponds to i?* , and that the energy density is proportional 
to fl gw cx (i?*i2*) 2 , with if* the Hubble time at the PT (c.f. Eq. © and Eq. ©). 

Afterwards, Hogan 1986 [l7| re-examines the problem in more detail, generalising the analysis to the EWPT as 
well, and providing for the first time an estimate of the shape of the GW spectrum. It is argued that the spectrum of 
GW energy density per logarithmic frequency interval must grow as fc 3 on frequencies lower than the characteristic 
frequency fc*, because the source is uncorrelated on time-scales larger than 1//3 and on length-scales larger than f?* 
(see section III Al) . From arguments based on the dynamics of the pressure waves generating the GW signal, [l?} also 
argues that the GW spectrum at frequencies higher than fc* should fall off at least as fast as 1/fc. 

A few years later, Turner and Wilczek 1990 [18[ study GW production by bubble collision in the context of extended 
inflation, where the transition between the inflationary and the radiation-dominated phases occurs through bubble 
nucleation and percolation. Again, it is found that the wavelength of the GW signal corresponds to the typical size 
of the bubbles when they collide, ii* ~ 1/if* (which in this case is of the order of the Hubble time at the end of 
inflation), and that the GW energy density is fl gw ~ 10~ 5 (f?*fi») 2 (c.f. Eq. ©, where p s = p to t in that case). 

In the above mentioned works, the GW signal was estimated using dimensional arguments and the quadrupole 
approximation, which gives the correct scaling but in principle breaks down for highly relativistic processes, for which 



8 



higher multipoles are also important. Kosowsky et al 1992 [lj| for the first time go beyond an analytical estimate of the 
GW signal, computing numerically the GW radiation produced in the collision of two broken phase bubbles in vacuum 
and Minkowski space. The result still depends only on the grossest features of the PT: the duration 1/(3, which here 
corresponds to the initial bubble separation and which is inserted as a time cutoff to break the 0(2, 1) symmetry of 
the problem; and the fraction of false vacuum energy liberated by the collision, i.e. the energy available in the source 
for GW generation (which in [IH - HH was taken to be of the order of the total energy density in the universe) . The GW 
energy estimated from the collision simulation (before red-shifting to today) is E gw /E vac ~ 10 _3 (-H*//3) 2 , confirming 
the usual dependence on the time/length-scale of the PT (c.f. Eq. ([S]), and note that this is also comparable to the 
factor (R*H*) 2 since in this case v\, = 1). The GW power spectrum shows a low frequency tail which is fitted in the 
paper to fc 2 8 (we know by causality that it should vary as k 3 at very small frequencies, see section Hi Ap . and an high 
frequency tail which is fitted as k~ 18 . The spectrum peaks at k pea k ~ 3.8 /3. Moreover, realises for the first time 
that, since it is the bulk motion of the bubble walls that produce most of the GW signal, it is possible to infer the 
GW emission by simulating only the motion of the bubble walls as a simple propagation of spherical, infinitely thin 
shapes instead of using the Klein-Gordon equation to evolve the scalar field. In this approximation, the gravitational 
radiation produced by bubble collision is due only to the TT part of the energy-momentum tensor of the uncollidcd 
envelope of the spherical bubbles, ignoring the interaction region. This 'envelope approximation' greatly simplifies 
the numerical simulation since it dispenses with the detailed dynamics of the scalar field and reduces the required 
computational power. It therefore allows to simulate the nucleation and collision of many bubbles, a more realistic 
model of the PT. 

This is first exploited in Kosowsky and Turner 1993 20], where the PT is simulated with the collision of 20- 
200 bubbles in vacuum. The nucleation rate is taken to be exponential in time with characteristic inverse time /?, 
r oc exp(/3i) (c.f. section Hi D[) . and the resulting GW spectrum agrees well with the two bubble case, except for a 
few features. First, the efficiency of GW generation is higher, leading to E gw /E vac ~ 0.06(iJ»//3) 2 . Second, the peak 
of the spectrum is at somewhat lower frequency, k pea k ~ 1.6/3, which corresponds to the peak of the energy weighted 
distribution of bubble sizes. Third, the spectrum has a slower high frequency decay, due to the late nucleation of 
smaller bubbles which is not present in the two-bubble case. However, the low numerical accuracy does not allow a 
conclusive statement on the slope of the GW spectrum at high frequency. 

Kamionkowski et al 1994 [21j generalises the analysis of [2fJ] to a PT happening in a thermal environment. The 
bubbles are described by means of the theory of relativistic spherical combustion, following the work of Steinhardt 
1982 [22]. The velocity of the bubble wall is therefore not any longer equal to the speed of light, and the energy 
momentum tensor sourcing the GW is no longer the one of the scalar field, but the one of the relativistic fluid: 
Tij — (p + p)viVj , where p and p are the energy density and the pressure of the relativistic fluid, and Vi is the velocity 
field describing the growth of the bubble. In Ref. 22] it is demonstrated that, in the case of bubble propagation by 
Jouguet detonations, is it possible to establish a very simple relation between the bubble wall velocity Vb and the ratio 
a of the vacuum energy density to the radiation energy density: 



/ \ VV3 + \Jol 2 + 2a/3 p vac 
Vb(at) = - , a=-r— , (10) 

1 + a Prad 

(where p* ad denotes the radiation energy density at the moment of the PT) . This formula is re-derived in [21( and it 
will be also widely used in the following literature on GW from PTs, which therefore implicitly assumes that bubble 
propagation is well described by a Jouguet detonation; however, this need not be the case, see section Hi Dl 

A second relevant quantity for GW production is first identified and evaluated in [2lj : the fraction k of vacuum 
energy that goes into bulk kinetic energy of the fluid as opposed to thermal energy. When the PT happens in vacuum, 
the liberated vacuum energy is transferred entirely into the propagation of the bubble wall (the gradient of the scalar 
field); for a thermal environment instead, the vacuum energy initially available is distributed during the PT among 
thermal energy of the plasma, bulk fluid motions at the bubble walls, and gradient energy of the scalar field. The 
latter need not be the dominant one, as in the case of vacuum bubbles; for moderately strong PT, it can actually be 
negligible with respect to the other components. The parameter k represents, in [2l| . an efficiency factor quantifying 
the amount of bulk kinetic energy which is the relevant source of the GW for a PT happening in a thermal bath. Using 
the velocity profile of a Jouguet detonation front, Ref. plj evaluates the energy momentum tensor TJy representing 
the propagation of one bubble, and provides a simple analytical fit of the efficiency factor as a function of a: 

. p Un 0.715a + 4/27 VW2 

K(a) — = - . (11) 

V ' Pvac 1 + 0.715a v ' 

Moreover, [21] analyses the spatial distribution of the energy momentum tensor and finds that for values of a > 0.01 
the stress energy is concentrated in a thin shell nearby the bubble wall. It is then concluded that, even in the case of 
a PT in a thermal bath proceeding through detonation, one can use the envelope approximation to simulate the PT 
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and evaluate the GW signal. This is how 21] proceeds to get the following formulas for the amplitude of the radiated 
GW energy density at the peak of the spectrum, and for the peak frequency {k peak denotes the physical wavenumber 
at the moment of the PT, while / = (fc/27r)(a*/ao) is the redshifted frequency today): 

^ gw \ ^10- 6 ^f^Vf^Vf V! aV-V, (12) 
9 l P eak \ (3 J \l + aj \0.24 + v$) V 9* J 

k peak * 1.9/3 , f peak * 5.2 • 10-3 mHz (A) (^) (^) * . (13) 

In this work as well, no conclusive statement is made on the slope of the GW spectrum at high frequency. 

The above formulas, together with Eq. (fTU)) and Eq. (|lip have been widely used in the following literature on the 
subject, because they provide a very simple estimate of the GW signal, which is completely determined once the 
parameters a, (3 and T* are known. However, the situation can be more complicated, since this simple result strictly 
applies only in the case of Jouguet detonations, and in the case in which the envelope approximation is well justified. 
For example, it does not account for friction acting on the bubble wall, which could slow down the wall propagation 
into a deflagration [23]. Assuming Jouguet detonations amounts to setting the value of the fluid velocity at the inner 
boundary of the velocity profile in the wall frame to the speed of sound. There is no justification for this choice, 
although it has been used in the literature for simplicity. Setting arbitrarily some boundary condition for the velocity 
leads to inconsistencies and corresponds to ignoring the constraints from the equation of motion of the Higgs field, in 



other words, corresponds to ignoring the friction term that restrains the bubble expansion. In Espinosa et al 2010 23 1 
it is demonstrated that k and Vb depend not only on a, but also on the friction 77, and numerical solutions for these 
parameters are provided, as well as analytical fits (in its appendix A). Moreover, the analysis of [23| allows to extend 
the definition of the parameter k to represent the full energy balance of the PT, including not only the bulk kinetic 
energy of the fluid, but also the gradient energy of the Higgs field (see section Hi D 1[) . We will see that this can be 
exploited to evaluate the GW signal continuously for any strength of the PT. In the following, we will adopt the 
results of [23[ for our analysis of the G W spectra in section III D 21 

An attempt to go beyond the envelope approximation and find a result applicable also to the case of deflagrations, 
where the bulk fluid motions are expected to be extended over a wider region around the bubble wall, has been 
carried on in Caprini et al. 2008 [24]]. Here, an analytical method is developed to determine the GW spectrum, which 
exploits the intrinsic randomness of the nucleation process: the GW spectrum is found by solving analytically Eq. ^) 
where the power spectrum of the GW source, given by the anisotropic stresses of the bulk fluid motions, is evaluated 
by Fourier transforming the two-point spatial correlation function of the velocity field describing the bubble. GW 
production from bubble collision is modeled in this approach by u sing Wick's theorem on the four-point velocity 
correlator which determines the anisotropic stress power spectrum [25[. The resulting GW power spectrum has a 
well defined shape. It is here re-established that at low frequency the spectrum must increase as k 3 , because the 
fluid motions are uncorrelated at large scales. The GW spectrum is found to peak at a wavenumber corresponding 
to the characteristic length-scale of the process, the average size of the bubbles and to decay at frequencies 
higher than the peak as k~ ls . However, in this work the anisotropic stress grows with time (as it should, because of 
bubble collisions), but it is abruptly switched off at the end of the PT. Caprini et al. 2009 [25| have shown that the 
peak and the slope at high frequency of the GW spectrum strongly depend on the time structure of the source, and 
inserting a discontinuity in time can lead to spurious effects which influence the shape of the GW spectrum. In [25| 
it is demonstrated that the correct time modeling of the GW source when accounting for bubble collision is that the 
source must be always correlated in time (completely coherent), and it must be switched on and off in a continuous 
but not diffcrcntiablc way. This leads to a peak in the power spectrum corresponding to the characteristic time of 
the source, i.e. the duration of the PT 1/(3 (as found in numerical simulations of bubble collision - for more details, 

c.f. HH). 

Shortly after the analytical estimation of [24T |. Huber and Konstandin 2008 [26[ simulated again the collision of 
bubbles in the envelope approximation, with considerably improved numerical accuracy, which allows to determine a 
larger portion of the spectrum and consequently a more careful analysis of the high frequency behavior with respect 
to [201 ] . The bubble wall velocity is treated as a free parameter, so the analysis does not rely on the assumption 
of propagation by Jouguet detonation (u& must only be sufficiently close to the speed of light for the envelope 
approximation to hold). The GW spectrum is found to increase as k 2 S at frequencies smaller than the peak (A; 3 at 
very small frequencies) , and to peak at a frequency which is somewhat smaller than the characteristic time-scale of the 
PT 1/(3, and slightly depends on the bubble wall velocity. The dependence of the peak amplitude with the parameters 
of the PT is broadly in agreement with Eq. (TT2|) . The most interesting result of [26| is that the GW spectrum is found 
to decay as 1/k at high frequency, considerably slower than previously found. It is argued that this behavior is most 
probably related to the many small bubbles nucleated at a later stage of the PT, and it is potentially very interesting 
for detection prospects. In section Hi D I we use the results of [26| to evaluate the GW signal by bubble collision. 
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2. MHD Turbulence 



GW production from turbulent motions has been first analysed in Kamionkowski et al 1994 [211 ] . where it is pointed 
out that bubble collisions cause a stirring of the primordial fluid on a characteristic scale given by the bubble size 
at collision leading to the formation of turbulent fluid motions with a Kolmogorov spectrum. The GW signal 
has been estimated in this work using the quadrupole approximation; the peak frequency corresponds to the inverse 
largest eddied turnover time tr ~ R^/vr, where vr is the typical velocity of the turbulent motions on scale R*. The 
high frequency slope is found to be fc -9 / 2 . The large-scale part of the spectrum is not calculated. 



The problem has been re-analysed later on in Kosowsky et al 2002 |27j and Dolgov et al 2002 [28(, which both 
go beyond the quadrupole approximation directly integrating the GW equation (|4]), and both remark that since the 
turbulence arises in a fully ionised plasma, it has to be treated as MHD turbulence: therefore one must add to the 
evaluation of the GW spectrum also the contribution from the turbulent magnetic field, besides the one from the 
turbulent velocity field. The procedure for the evaluation of the GW spectrum in 27] is fully analytical (as a matter of 
fact, there have been no numerical simulations of GW production by MHD turbulence to date, contrary to the bubble 
collision case). Dolgov et al. 2002 [28| uses the result of [I?) for the GW spectrum, but parametrically introduces also 
power-law deviations from the Kolmogorov spectrum, and generalises the analysis to the case when the time during 
which turbulence is active is of the order of the Hubble time. Both analyses confirm the finding of [2l[ in what concerns 
the peak frequency, k pea k — 1/tr, but the high frequency slope that they find is different: fc -7 / 2 in the Kolmogorov 
case. However, both these analyses have some fundamental problems. First of all, they only evaluate the GW signal 
arising from the Kolmogorov part of the turbulent velocity spectrum, which is that part of the spectrum with slope 
fc— 11/3 m a range of wavenumbers comprised between the stirring scale i?„ and the dissipation scale. However, the 
power spectrum of the velocity field does not vanish on scales larger than the correlation scale (corresponding to the 
largest eddy size i?*), so the turbulent, as well as the GW spectrum, must cover the entire frequency range. Moreover, 
the extension of the Kolmogorov slope fc" 11 / 3 up to the stirring scale i?* causes an overestimation of the turbulent 
kinetic energy density, and therefore of the GW signal at the peak frequency. The second problem is that in these 
references the source is modeled as stationary in time, in a way that the two point correlator of the source at different 
times only depends on the time difference. However, by definition the MHD turbulence must decay in time, because 
the source of stirring (bubble collisions) is not eternal. Consequently, in this cosmological context, it is not fully 
justified to assume a stationary source. The third and most important problem, is an error in the dispersion relation 
of the emitted GW, which instead of satisfying the usual dispersion relation, by an artefact of the calculation ends 
up having the same dispertion relation of the turbulent source, i.e. u) cx fc~ 2 / 3 . The formulas for the characteristic 
peak frequency and the GW amplitude at the peak based on these analyses have been collected in Nicolis 2004 29] , 
and have been widely used in the following literature: 

^U-io- 4 4.^ 2 (f) 2 (^)\ (14) 

W * ^ , f P ea k * 3.4 • 10-VHz ^ ( A )( -^-) (±.) 1 . (15) 

p it,' p v b \hJ viooGevy VlOO/ v ' 

An attempt to correct the above mentioned problems has been carried on in Caprini and Durrer 2006 t 30] . Here it is 
demonstrated that the large scale part of the turbulent velocity spectrum is indeed present, and it is characterised by 
a Batchelor spectrum (i.e. it grows as fc 2 ): this allows to evaluate the full GW spectrum, not only the high frequency 
slope. The time dependence of the turbulent source is modeled not as stationary but based on an heuristic model of 
the oscillations of the random velocity field. This model is inspired by Richardsons energy cascade through eddies 
of different sizes, motivated by the fact that the time dependence of the turbulent flow is not lost in the ensemble 
averaging. Moreover, the correct dispersion relation for the free propagation of the GW is assured, i.e. uj = k after 
the turbulent source ceases to be active (for more details, see Caprini et al. 2006 [31]). The GW spectrum found 
in this work is quite different from the result of [27L [28| , with a peak frequency that corresponds to the size of the 
bubbles at collision k pea k — 1/-R*, a low frequency slope as k 3 by causality, and a high frequency slope as fc~ 8 / 3 . 
The GW spectrum from the magnetic field has the same characteristics but a different high frequency slope, fc~ 5 / 2 , 
due to the different assumption for the time evolution of the magnetic and turbulent sources. The analysis ends 
with a comparison of the efficiency of the turbulent and magnetic sources, finding that since the magnetic field is 
a coherent source and it is active for a long time, there is a very effective conversion of magnetic energy in GW 
energy at horizon crossing. However, in this reference as well the Kolmogorov slope is extended up to the largest 
eddy size, so the amplitude of the GW signal at the peak is overestimated. This problem persists also in the analysis 
of Gogoberidze et al. 2007 [32j, which keeps the assumption of stationarity of the source, but corrects the problem 
of the dispersion relation. The results of [30| have been used in Megevand 2008 [H[ to evaluate the GW signal in 
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the case of deflagrations, and it is shown that for some values of a < 1 and of the PT temperature, one can reach 

^gw 00 10 

The analysis of Kahniashvili et al. 2008 [34, 35] shows that helical turbulence, which could be produced during a 
first-order PT if there is macroscopic violation of parity (37|, generates circularly polarized GW. Kahniashvili et al. 
2005 ;36( concludes that it is unlikely that this GW signal has sufficient amplitude to be detected by the planned 
GW detectors at that time, but this conclusion is changed by the subsequent analysis of [34], [35[, which also accounts 
for the presence of magnetic fields and for the time evolution of the MHD turbulence by inverse cascade. The MHD 
turbulence is assumed to be stationary and decorrelating in time, but an overall time dependence is inserted through 
the decay of the MHD spectrum with the turbulent cascades. It is found that, compared to the non-helical case, the 
spectrum of GW generated by MHD helical turbulence peaks at lower frequency with larger amplitude and could have 
been detected by LISA. However, here as well the Kolmogorov slope is extended up to the largest eddy size, leading 
to an overestimation of the signal. In Caprini et al. 2009 [38j . the same analysis has been carried on with similar 
models of the inverse cascade, but in the MHD turbulence spectrum the transition from the low wavenumber part of 
the spectrum to the inertial range is modeled in a more realistic way, as opposed to extending the Kolmogorov slope. 
This leads to an amplitude of the GW signal which is lower by nearly one order of magnitude. In Ref. [38| it is argued 
that the different assumptions for the unequal time correlation of the source may also play a role in determining this 
difference (in [34|, HH the source is assumed to be stationary with exponential decorrelation, while in [38[ the source 
is assumed to be completely coherent, which leads to an underestimation of the GW signal). 

The last analysis of GW emission from MHD turbulence generated during a first-order PT and freely decaying 
afterwards is the one of Caprini et al. 2009 [nj]. This paper takes into account only non-helical MHD turbulence, but 
models the source in a time-continuous fashion, starting with zero kinetic energy and building up the Kolmogorov 
spectrum after one eddy turnover time (instead of instantaneously), and then accounting for the free decay of the 
source (since the stirring due to the PT lasts only for about one eddy turnover time). The source power spectrum 
is again modeled with an interpolation in wavenumber between the large and small scale behaviours, determined 
respectively by causality and by the Kolmogorov (or Iroshnikov-Kraichnan) theory. The time de-correlation of MHD 
turbulence mimics the model proposed by Kraichnan, as in [34, 35]. The claim of this analysis is that, in the case of 
MHD turbulence, neither the coherent [38| nor the stationary [34 , 1 3 51 ] approximations previously used in the literature 
are the correct ones, but the time dependence of the source has to be modeled using a top hat ansatz, in order to 
best reproduce the Kraichnan decorrelation model. Moreover, it is demonstrated that it takes several Hubble times 
for the MHD turbulence to decay, due to the very low viscosity: therefore, MHD turbulence has to be modeled as a 
long lasting GW source, which causes an amplification of the GW spectrum at very large scales by about two orders 
of magnitude with respect to the short lasting case. The GW spectrum found in this work shows a peak frequency 
that corresponds to the size of the bubbles at collision k pea k — 1/i?*, a low frequency slope as fc 3 by causality, and a 
high frequency slope as fc~ 5 / 3 if the source has a Kolmogorov spectrum, and as fc _3//2 if the source has an Iroshnikov- 
Kraichnan spectrum. The analytical evaluation maintains a certain level of intrinsic uncertainty which probably can 
only be addressed by numerical simulations of relativistic MHD turbulence. We nevertheless use the result of [Tl| in 
our estimation of the GW signal in the eLISA frequency band, section Hi D[ 

C. Overview of the Literature: GW Production in Models of the EWPT 

Parallel to the analyses mentioned above, which aim to find the form of the GW spectrum by bubble collision and 
MHD turbulence in terms of a few free parameters describing the source, there have been a number of works which 
consider specific models of the EWPT and evaluate the expected GW signal in the context of the given model. This 
allows to determine, for example, the parameters a and /3/-H* in connection with the symmetry breaking potential, 
and from this make definite predictions of the GW signal. Often the formulas of Eq. (fT2"|) and Eq. (|14l) have been used, 
more or less implicitly assuming propagation by Jouguet detonation, the envelope approximation, the validity of the 
Kolmogorov spectrum up to the relevant eddy scale, and not correctly accounting for the GW dispersion relation in 
the case of turbulence. 

Apreda et al. 2001 HH and Apreda et al. 2002 [U use Eq. (TJ2]) and Eq. (0, together with Eq. (rU]) and 
Eq. (|11[) . to evaluate the GW production during the EWPT in different supersymmetric models where, contrary 
to the case of the standard model (SM), the transition can be first-order. The GW energy density is found to be 
low in the minimal supersymmetric standard model (MSSM), but in some parameter range of the next to minimal 
supersymmetric standard model (NMSSM) can be as high as £l gw ~ 4 • 10~ . This happens since the cubic terms 
present in the NMSSM tree-level Higgs potential can lead to a much stronger PT. This analysis has been redone 
in Huber and Konstandin 2008 [4lJ for the nMSSM case (nearly minimal supersymmetric standard model), and 
considering also the SM augmented by dimension six operators (which can be viewed as the low energy effective 
description of some strongly coupled dynamics at the TeV scale). In this model, the barrier in the Higgs potential 
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is also present at tree-level [42f. The method used to determine the bubble configurations is improved with respect 
to [4(| (where it led to an overestimation of the signal). In [4lJ, Eq. (fT2l is used for the estimation of the GW 
signal by bubble collision, while for MHD turbulence both the results of [28l| and those of [3(| are used. Eqs. (fT0|) 
and Eq. (fTTj) make the connection between the bubble wall velocity and the strength of the PT (implicitly assuming 
Jouguet detonations). The strong point of this work is that the parameters a and /3/H* are evaluated directly from 
the three dimensional Euclidean action S3 of the symmetry breaking theory, and are therefore directly related to 
the features of the potential. The connection is such that large values of a imply small values of f3/H*, because 
stronger PTs happen at lower temperature and create larger bubbles. This is bad for detection, since the peak of 
the spectrum is shifted to lower frequencies, out of the sensitivity range of LISA (see also the discussion in section 
III Dp . As a consequence, [4l[ finds that the GW signal is in a large portion of the parameter space too small to be 
detected. This result is in qualitative agreement with Delaunay et al. 2008 [Hj], which provides a more detailed and 
complete one-loop finite temperature calculation of the Higgs potential augmented by dimension-6 operators, as well 
as the corresponding GW signal. 

Another model leading to a strong first-order PT is provided by extensions of the SM with hidden sector scalars 
coupled to the Higgs (Espinosa and Quiros 2007 Espinosa et al. 2008 [Hj]). The GW signal from the PT in the 
context of this model is calculated, by means of Eq. (fT2"|) . Eq. (ITU1) and Eq. (fTTj) . but it is found that no observable 
traces are expected for LISA. Kehayias and Profumo 2010 (46[, on the other hand, calculates the GW signal from 
a generic SM-like potential: using a semi-analytic approximation for the Euclidean three-action S3, approximated 
expressions are derived for the tunneling temperature, a and fi/H*, as functions of the parameters in the potential, 
which allow to fully specify the GW spectrum. Singlet, triplet, higher dimensional operators and top-flavor extensions 
to the Higgs potential are considered, and it is found that the addition of a temperature independent cubic term in the 
potential, arising from a gauge singlet, is the most promising model to enhance the GW signal. The main difference 
with the previous analyses is the application of the semi-analytic approximation for the three-action. 

More promising for detection are GW arising from PTs in the context of nearly conformal dynamics at the TeV 
scale, and, by holography, of extra-dimensional warped geometries. This has been first proposed in Randall and 
Servant 2007 (47] . where it is demonstrated that the Randall-Sundrum 1 (RSI) model can provide a strong signature 
in GW, generated during the first-order PT from the AdS-Schwarzshild phase to the RSI geometry, which proceeds 
through the nucleation of brane bubbles. In this work as well, the GW signal is estimated using Eqs. (fT2l[T4)) and (|T0l 
Hip , reducing the problem to the unique determination of the parameters a and /3/i?*. It is found that the PT can be 
very strongly first-order and gives an interesting GW signal in the LISA frequency band. Konstandin et al. 2010 (48j 
confirms this finding, improving the analysis on the model side, by partially accounting for the back-reaction effect 
from the scalar field undergoing the transition on the metric, and on the GW estimation side, using the formulas given 
in the most recent simulations of (26|. In Konstandin and Servant 2011 [4!| it is emphasised that a detectable GW 
signal would be a smoking gun of nearly conformal dynamics at the TeV scale. The reason is twofold: on one side, 
due to the large supercooling, these models give rise to GW signals which are typically larger in amplitude than for an 
EWPT arising from a generic polynomial potential (for which a fine-tuning of parameters is required to have a large 
supercooling); on the other side, the PT results from dynamics at the TeV scale, and therefore the peak frequency is 
about one order of magnitude larger than for a generic EWPT, which improves the prospects for observation. 

We conclude our review of the literature with the work of Nicolis 2004 [29| and Grojean and Servant 2007 (5C| : 
these papers aim to provide model independent analyses of the GW signal from a first-order PT, scanning the values 
in the parameter space (a, /3/H*) which lead to an observable signal. However, these analyses need to be updated in 
at least two aspects. First of all, they use Eqs. (IT2"llT4"f and (ITUllTTj) to estimate the GW signal: as we have mentioned, 
these equations are only valid in the case of Jouguet detonation and in the envelope approximation, and give a wrong 
estimate of the GW signal from turbulence. Secondly, they consider a and /3/-H* as independent parameters, while 
it is demonstrated, for example in Ref. [4l|, that this is not the case in the context of complete EWPT models. In 
the following, we aim to provide the most reliable possible analysis of the GW signal relevant for eLISA, given the 
present status of the research on this subject, correcting for these aspects which are biasing previous analyses. 

D. Possible GW Signal in the eLISA Frequency Band 

The GW signal from a first-order EWPT is given by the sum of the GW spectra due to bubble collision and MHD 
turbulence. Concerning bubble collision, Huber and Konstadin 2008 (2(| provides the most recent simulation of the 
collision of bubbles in the envelope approximation. We therefore take the result of this work for our estimation of the 
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signal in the eLISA frequency band. The GW spectrum resulting from the simulation is [26j, |48| 

h»Q ml - 167- IP' 5 (^Y(J2-Y( ° A1V " ) f^Y 3.8/^/" 



This result is comparable to the one of Eqs. (fT2"j) and (|13p . with the difference that the peak frequency depends 
somehow on the bubble wall velocity v^, and also that the dependence of the GW amplitude on Vb is slightly modified. 
The peak frequency approaches the one of Eq. (|13l) in the limit v — > 0, while it is slightly smaller for large velocities 4 . 
This corresponds to a peak wave- number of the order of (3, the inverse characteristic time-scale of the PT. As mentioned 
before, the physical reason is that the bubble collision source can be approximated as totally coherent (see [TBI. [25j). 
The most notable aspect of the analysis of (2(| is that, thanks to the improved numerical accuracy, it has been possible 
to determine the shape of the GW spectrum. In Eq. (fT6|) we have introduced the fit to the shape of the GW spectrum 
resulting from the numerical simulations, given in Eq. (19) of (26[. Note that this fit is optimized for a frequency range 
close to the peak frequency, and does not correctly reproduce the asymptotic low frequency behavior, which should be 
/ 3 by causality (see section liTA"f . We use it nevertheless to determine the full spectrum, since anyway the interesting 
frequency region for detection is the one close to the peak. The above result is also in broad agreement with the 
rough estimates given in Eqs. © and ©. In particular we see, from ©, that for the bubble collision case the relative 
energy density available in the source for the GW generation is given by {p 8 / ptotj* = Pkin/ '{p*ad~ >r Pvac) = Ka/(l + a). 

Concerning MHD turbulence, for our estimation of the signal in the eLISA frequency band we use the result of the 
most recent analysis of Caprini et al. 2009 ITll . Here the GW spectrum is evaluated analytically; a simple fit to the 
result of [Tl| is given in Caprini et al. 2010 [14J and reads: 
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where fc» = 27r / a a /a* is the physical wave-number at the time of production. Since these formulas are less commonly 
used in the literature than Eqs. (TT6"|) and (fTT)) . let us analyse them in some detail. In [ll| it is argued that for MHD 
turbulence, both the velocity field and the magnetic field contribute to the GW source: the collision of the bubbles stirs 
the fluid leading to kinetic turbulence, and the kinetic energy of the turbulence is then converted partly into magnetic 
energy by the MHD cascade, until equipartition is reached, i.e. the same amount of energy is contained into turbulent 
motions and the magnetic field, px — Pb- Because of the (well justified) assumption of equipartition 5 , we refer only to 
the kinetic energy density of the turbulent source to represent the total energy available in the MHD source, setting 
Pt + Pb — %Pt- In [3, the relative energy available in the source for the GW generation is given in terms of the 
parameter 51* = (p s / Ptot)* — ^Pt/p*', in Eq. (fTS)) . we have expressed this in terms of Ka/(l + a) = Pkin/(p* a d + Pvac)- 
As for the bubble collision case, we adopt the parameter k quantifying the kinetic energy of the fluid to represent 
the total energy available in the MHD source. We are therefore assuming that the efficiency of converting the bulk 
kinetic energy of the fluid motions due to the bubble expansion into turbulent motions after collision is just one: all 
the energy of the bulk motions pkin is converted to turbulent kinetic energy and is then transformed partially into 
magnetic energy by the cascade, pkin — 2pr- A proper estimate of the efficiency with which the bulk motions due to 
bubble expansion transform into chaotic turbulent motions has never been performed; however, the kinetic Reynolds 
number of the universe at the time of the EWPT and at the scale of the bubble radius is of the order of Re(i?„) ~ 10 13 
[ll| . Given this extremely high value of the Reynolds number, corresponding to an extremely low value of the kinetic 
viscosity, it seems very reasonable to consider that the bulk fluid motions are converted into kinetic turbulence without 



any dissipative loss, at least at the scale of the bubble radius. Furthermore, in Ref. [11| it is demonstrated that, also 



4 In principle, in the limit of small bubble wall velocity the envelope approximation breaks down, because one enters the regime of 
propagation by deflagration. 

5 The tail of the MHD spectrum is not firmly established in the theory of MHD turbulence, so here we consider for simplicity that the 
magnetic field as well develops a Kolmogorov tail in its spectrum, just as the kinetic turbulence: fc — 11 / 3 . In the Iroshnikov-Kraichnan 
shape was also considered, but the difference is irrelevant as far as detection prospects are concerned. 
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due to the very small kinetic viscosity, the MHD turbulence dissipates slowly in the early universe: therefore, the 
GW source is in general active for many Hubble times after the completion of the PT. This fact needs to be taken 
into account when evaluating the GW spectrum. That is why, the relative energy density available in the source for 
the GW generation has been normalised in fill ] to the total radiation energy density in the universe, f2* = (p s /ptot)* 
with ptot = Prad after the completion of the PT. Note that this is not in contradiction with what we did above, since 
after the completion of the PT ptot — Prad, while during the PT ptot = p* aa - + Pvac (we remind that p* ad denotes the 
thermal energy density present in the universe at the moment of the PT) . For the pre- factor in Eq. (|18j) we are simply 
accounting for the fact that the vacuum energy plus thermal energy present at the moment of the PT is completely 
converted into radiation energy at the end of the PT (plus, obviously, the turbulent kinetic and magnetic energies, 
but these are subdominant with respect to the radiation energy density in order not to spoil the homogeneity and 
isotropy of the universe) . 

The fact that the source is long lasting changes in Eq. (|18p the dependence of the GW energy density with (p s /ptot)* 
and H*//3 with respect to the scaling given in Eq. (O, and causes the slope of the MHD signal to turn from / 3 to f 2 
at sub-horizon scales a*iJ*/(87ra ) < / < f pea k [llj . Moreover, in the case of MHD turbulence, the peak frequency 
corresponds to the inverse characteristic length-scale of the source, i.e. the bubble size towards the end of the PT, 
i?* = Vb//3. Again, the physical reason for this is related to the properties of temporal correlation of the GW source: 
in [ll[ is it argued that MHD turbulence should be modeled with a top hat temporal correlation, since this is the best 
way to reproduce the Kraichnan decorrelation model, which state that the turbulent velocity field on a given scale I 
is not any longer correlated after a characteristic time-scale given by the eddy turnover time on the scale £, Tg — vgj I. 
GW sources with a top-hat decorrelation structure typically lead to GW spectra peaked at the inverse characteristic 
length-scale of the source [HI, [25| . 

From Eqs. (IT5|) . (fTT)) and (fT5)l and P^|) it appears that five parameters must be known to determine the GW 
spectrum, namely a, j3/H*, k, v and T». These are not independent parameters, and their evaluation, once a specific 
model for the first-order PT is given, requires involved numerical computations: a, fi/H* and T* are in principle 
determined from the shape of the potential at finite temperature, while for k and Vb one needs to know the details of 
the propagation of the bubble wall in the surrounding plasma. We now revise how the problem of evaluating these 
parameters has been tackled in some works in the literature, and explain the strategy we adopt to use the results of 
these analyses to calculate reliably the GW spectra. 



1. Phase Transition Parameters 



In the first a naly ses of the GW signal from PTs, a general approach was provided to calculate the parameters /3/H* 
and T* [l9l.l20l.l52f. The approach is based on the Taylor expansion of the tunneling action around the time when the 
transition completes: 



S(t) = S(t*)-p{t-t*) + G{(t-U) 2 ), so that -J- = T*% 

H* dl 



(20) 

T, 



with dT/dt = —TH. The 'temperature of the transition' T» is defined as the temperature at which the probability of 
nucleating one bubble per Hubble volume per Hubble time is equal to one: denoting T(t) the bubble nucleation rate 
per unit volume and time, one has 

r(f) ~ T 4 e" s « , -> S(T*) = -41n^, (21) 

mp 

where mp is the Planck mass. Since the Hubble time is the characteristic time-scale, Refs. set (3 ~ H*S* as an 

order-of- magnitude estimate, so that /3/H* ~ —4 ln(T*/mp). This provides a relation among the parameters j3 and T*, 
and for the EW scale gives the value P/H* ~ O(10 2 ). Using this relation reduces the number of required parameters 
to determine the GW spectrum from five to four; furthermore, under the assumption of Jouguet detonations, this 
can be combined with Eqs. (ITU1) and ([TT|) . which give the bubble wall velocity and the efficiency factor as functions 
uniquely of a. Therefore, in this approach the parameter space to fully specify the GW signal can in principle be 
reduced to only two parameters: either the couple (a, T*) or the couple (a, /3/H*). As already mentioned, this 
has been exploited for example in Nicolis 2004 (29| and Grojean and Servant 2007 [50]. Since then, there have been 
other works which have analysed the problem using more accurate estimations of /3/H*. Note, however, that this 
necessarily implies a loss of generality, because one has to focus on a specific model of the PT. 

A more accurate way to proceed given a model, i.e. some extension of the SM of which one wants to study the 
EWPT, is explained for example in Huber and Konstandin 2008 [4l[ and Espinosa et al. 2008 [45|. In general, 
one needs to determine the bounce solution of the three-dimensional Euclidean action S^iT), which quantifies the 
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probability of tunneling through T ~ T 4 [S 3 (T)/2ttT] 3 / 2 exp(-S 3 (T)/T). One can then calculate the average number 
of bubble nucleations per Hubble volume P(T), and the fraction of space that is covered by bubbles neglecting overlap, 
f{T): 

w-r?®, /(n^rtf^f), (22) 



T 5 4 ' v ' 3 J T T H 

where T c denotes the critical temperature of the PT and R the bubble radius. This allows one to determine precisely 
the initial and final temperatures of the PT: the temperature at which the PT starts is defined by the nucleation of 
the first bubble, P{T n ) = 1, which leads to the well known condition S 3 (T n )/T n ~ 142; the temperature at which the 
PT ends is defined as f(Tj) — 1. Moreover, knowing the three action as a function of temperature, one can calculate 
P/H* = T d{S^/T) / dT , which is in general a function of temperature and has to be evaluated towards the end of 
the PT to represent, as a matter of fact, the 'duration' of the PT [4l|. Alternatively, it is possible to relate /3 to the 
typical bubble size at the end of the PT through Vf,, (R) — 3vb/P(T), where (R) can be estimated from the maximum 
of the bubble volume distribution and has to be evaluated at an intermediate temperature between T n and Tf [4l| . 
Without going into many details, here we only want to point out that from the point of view of estimating the GW 
signal, the importance of adopting the more accurate way to evaluate /3/H* described above resides in the fact that 
this reveals a general relation among the parameters /3/H* and a: large values of a, i.e. strong PTs, lead in general 
to small f3/H* and therefore to a lower peak frequency. This is due to the characteristic shape of the three action as 
a function of temperature (see e.g. Fig. 5 of [45]), and has relevant consequences as far as detection prospects are 
concerned: while large a make the amplitude of the signal increase, the fact that this corresponds to small /3/iZ* 
shifts the peak of the spectrum to lower frequencies, progressively out of the detectable frequency band (c.f. Fig. [5] 
and also the discussion in (4l|V 

We now turn to the analysis of the parameters Vb and k, which together with a and /3/-H* are the other two 
parameters that one would like to predict for a given particle physics model leading to a first-order PT at temperature 
T*. As already mentioned, Eqs. (fTUj) and (fTT|) do not hold in all generality, because depending on the characteristics 
of the PT the propagation of the bubble in the surrounding plasma need not be a Jouguet detonation. Other kinds of 
solutions for the bubble expansion are possible, such as weak detonations, deflagrations, hybrids and runway solutions. 
While in detonations and deflagrations the wall velocity is supersonic or respectively subsonic, leading to a rarefaction 
wave behind the bubble wall or a compression wave in front of it, in hybrid solutions both the rarefaction and the 
compression waves are present, and in runaway solutions the wall continuously accelerates, nearly without interaction 
with the surroun ding plasma. Model independent studies of the hydrodynamics of bubble growth are performed for 
example in Refs. |23l l5ll [53j . We concentrate on the work of Espinosa et al. 2010 (23[, where it is shown that the 
knowledge of the parameter a alone is not enough to determine the bubble propagation mode. For a given a, there 
are many possible values of the wall velocity Vj, and of the efficiency coefficient k: to fix the solution, one has to specify 
the friction acting on the bubble wall because of its interaction with the surrounding plasma. The friction grows with 
the velocity of the expanding bubble and opposes the driving force due to the difference of the free energy across the 
wall, until the bubble growth reaches a steady state at a constant terminal velocity. The hydrodynamical analysis 
alone cannot completely determine the terminal wall velocity v^. Rather than solving the complicated system of 
Boltzmann equations describing the distribution functions of all particle species and their interaction with the Higgs 
field, as done for example in [54. [55[ ■ Ref . [23| uses a phenomenological description of the friction 77 as a parameter 
independent of Vb (see also [53[~ 5(J~ El]] ) • Model-independent contour plots for Vb and k are then given in the (a , 77) 



plane: from a given particle physics model, one can in principle evaluate a and 77, and once these two quantities are 
known, the results of [23| provide a way to determine Vb and k. 

Moreover, Espinosa et al. 2010 [23] evaluates the energy budget of the PT, i.e. the efficiency of the transfer of the 
vacuum energy into gradient energy of the Higgs field, bulk kinetic energy of the fluid motions at the bubble wall, and 
thermal energy. This is extremely important, since it allows to estimate the relative energy available in the source 
for the GW generation realistically, given a value for the parameter a. In the following, we use the results of [23| for 
the first time in the calculation of the GW signal. As far as production by bubble collision is concerned, GW can be 
sourced by both the gradient energy of the Higgs field and the bulk kinetic energy of the fluid, since both components 
lead to tensor anisotropic stresses. The efficiency parameter k, appearing in Eq. (|16p . needs therefore to be divided 
into two components: k = + n v , where denotes the fraction of vacuum energy that goes into gradient energy 
of the Higgs field, while n v denotes the fraction that goes into bulk kinetic energy. Inserting the total k in Eq. (|16l) . 
one can use this equation to model in a continuous way the GW signal from both weak (for which n v is important) 
and very strong (for which is important) PTs. One has then + k v + pther/Pvac = 1, where pther is the thermal 
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energy liberated by the PT, and it is the only component that does not source GW 6 . On the other hand, the gradient 
energy of the Higgs field does not, by definition, source the MHD turbulence: the efficiency parameter n appearing 
in Eq. (|18|) corresponds only to k v . 

The parameters kq, k v and pther/pvac are calculated in [231 ] as a function of a and for two different values of the 
friction, 77 = 0.2 and 77 — 1 (c.f. Fig. 12 of [23]). As the PT gets stronger (a increases), the total kinetic energy 
Pkin + Pkin naturally gains importance with respect to the thermal energy pther- Ref- [23| finds that, as long as the 
bubble expansion proceeds as a deflagration or a detonation, the gradient energy in the Higgs field is completely 
subdominant and n = n v . But when a overcomes a certain value, the bubble propagation enters the runaway regime. 
This means that the PT is so strong, that the friction exerted on the bubble wall by the plasma is not sufficient to 
counterbalance the driving force due to the pressure difference in the broken and symmetric phases, and the wall 
continuously accelerates. When this happens, is no longer negligible. The value of a after which one enters the 
runaway regime depends on the friction, and corresponds to a > 0.25 for 77 = 0.2, and a > 1 for rj = 1: as expected, 
[23^ finds that as the friction increases, the ratio between k v and increases for fixed a, i.e. more energy is transferred 
to bulk motions of the fluid. On the other hand, for fixed friction and increasing a, the gradient energy of the Higgs 
becomes larger and lager with respect to the bulk kinetic energy. 

The analysis of Ref. |23[ leads to the conclusion that only for a restrained set of values for a and 77 the propagation 
solution is a detonation, while most likely solutions are deflagrations or runaway solutions. This needs to be taken 
into account when estimating the GW signal. Even though the value of the friction 77 is very difficult to determine in 
a given EWPT model (it is in fact known only in the SM [H[ and in the MSSM (58|), still the analysis of Ref. [23| 
provides a much more realistic estimate for the parameters Vb and k than what is usually assumed using Eqs. (| 10[) and 
([TT]) : we therefore use the result of this analysis to evaluate the GW signal. Moreover, introducing k^, + k v in Eq. (|16p . 
and k v in Eq. (fT5|) . we can add the bubble collision signal and the MHD turbulence signal in a fully consistent way 
for any given value of a. 

2. Estimate of the GW Signal 

From the discussion in the above section, it appears that one is left with four parameters that fully determine the 
GW signal, according to Eqs. (|16H19[) : a, /3/i2*, 77 and T*, since k and Vf, can be determined from a and 77 [23]. For 
a given model of the PT these parameters are all known: a, and T* from the finite temperature potential, 

77 from the particle content and the interactions of the theory. Rigorously, these parameters are not independent of 
each other, and one cannot in principle perform a completely model-independent study of the GW signal letting the 
parameters vary freely. Nevertheless, at the end of this section we present the detection prospects of the GW signal 
in the eLISA frequency band from PTs in the GeV-TeV energy range, trying to remain as model-independent as 
possible. We have also chosen two examples presented in the literature, which turn out to be the most promising for 
detection with eLISA and for which we calculate the GW spectra. 

First of all, let us consider the analysis of Huber and Konstadin 2008 [4l|. In this work, two models of the EWPT 
are studied: the nMSSM and the SM augmented with dimension 6 operators (c.f. section lTl B[) . In the latter model, the 
PT can be considerably stronger with respect to the nMSSM, leading to a higher GW signal: we therefore choose this 
model in the present analysis. Ref. [4l[ numerically evaluates the parameters of the PT a, i?*-ff* and T„, for different 
values of the parameters in the model: the suppression scale of the dimension 6 operator M, and the quartic coupling 
A - related to the physical Higgs mass. It then estimates the GW signal using the formulas presented in section Hi B I 
here we update this analysis by employing Eqs. p^]) - p^]l and the work of Espinosa et al 2010 [23] to relate k and Vb 
to a and rj. The results are given in Fig. [2] for two different values of 77. These figures are obtained as follows: each 
GW spectrum corresponds to a given set of parameters a, R^H* and T* (resulting from a choice of the energy scale 
M), according to table I of [4l|. To translate the parameter R^H* into P/H* one needs to know the bubble velocity 
Vb, since ~ 3vb/f3 [4l[: for each value of a, and for two fixed values of the friction 77 = 0.2 and 77 = 1, we obtain 
the wall velocity Vb(a, 77) by linear interpolation of the contour plot values given in Fig. 10 of Ref. [23]. This allows a 
quite general estimation of Vb which goes beyond the usual assumption of Jouguet detonations. Knowing Vb, we can 
relate = 3vb/ (i?*-H*) directly to a: again, we do a linear interpolation of the values of i?*-ff* given in table 1 of 

[23| for different values of a. We can therefore determine, for a couple of values (a, 77), the corresponding ft/H* which 



To summarise our notations, the total energy density in the universe through the PT can be expressed as ptot = P* = pvac + P ra d 
just before nucleation, which becomes p va c + P* a[ j = Pkin + Pkin + Pther + P* ra d during bubble nucleation, and which becomes then 
Pkin + Pkin + Pther + P ra d = Prad a t the end of the PT (after percolation, and neglecting the remaining energy density in MHD 
turbulence). Furthermore we have set a = p V ac/p* rad , «0 = pf. in /p V ac and k v = p v kin / pvac- 
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SM + dimension six operator, j]=0.2 SM + dimension six operator, 77= 1 




f (Hz) f (Hz) 

FIG. 2: GW spectra for the EWPT in the SM augmented by dimension six operators [Jl[ . Left panel, 77 = 0.2, Right panel, 
77 = 1. The red line shows the eLISA noise curve. The GW spectra correspond, from bottom to top, to increasing values of a 
as given in Table 1 of Huber and Konstandin 2008 0: a = 0.128 , a = 0.201 , a = 0.311 , a = 0.586 , a = 1.197 , a = 2.268. 
The values of /3/H* and T* decrease accordingly (see discussion in the main text). 



one expects in the EWPT model under analysis, i.e. SM model extension by dimension 6 operators. To know the 
GW signal, it remains to evaluate the efficiency factor for each value of a, and for the two values of the 

friction 77 = 0.2 and rj = 1, we determine it from figure 12 of [23]. As explained above, the GW signal from bubble 
collision, given in Eq. (|T6]) . is due to the total k = k,^ + n v , since both the gradient energy of the Higgs field and the 
kinetic energy of the bulk fluid motions lead to anisotropic stresses when the bubbles collide. On the other hand, 
only the kinetic energy of the bulk fluid motions contributes to the GW signal from MHD turbulence: in Eq. (IT51) we 
therefore substitute n — n v . As mentioned above, we assume that after bubble collision all the bulk kinetic energy is 
transferred to MHD turbulent energy without any dissipative loss, due to the extremely low kinetic viscosity of the 
plasma. 

Fig. [2] shows the sum of the GW signal from bubble collision and MHD turbulence for rj = 0.2 and 77 = 1. The 
peak of the GW spectrum from MHD turbulence occurs at higher frequencies, since it corresponds to the inverse 
characteristic length-scale of the PT, the bubble radius at collision i?* (Eq. (fT9|) . as opposed to Eq. (jTTJ) ) . It appears 
that as a increases, the signal from bubble collision gain s importance with respect to the one from MHD turbulence. 
This can be explained as follows (see also Fig. 12 of [23|]): in the deflagration or detonation regime, the same kinetic 
energy given by k v determines both signals, from bubble collision and from MHD turbulence. But as soon as one 
enters the runaway regime, the gradient energy of the Higgs field starts to play a role, and eventually, as the PT 
gets very strong, it completely dominates over k v . When this happens, there is no more GW production by MHD 
turbulence and the signal is only due to bubble collision. Clearly, the value of a for which this happens depends on 
the friction: if the friction is large, the transfer of the vacuum energy to bulk kinetic energy is more efficient and one 
needs a stronger PT (larger a) to reach the point in which MHD turbulence is not any longer present. Therefore, in 
Fig. [2] the GW signal from MHD turbulence is still visible at large a ~ 2.27 for 77 = 1; on the other hand, for 7/ = 0.2, 
when a > 0.5 the gradient energy of the Higgs field becomes dominant with respect to the bulk kinetic energy, and 
no MHD turbulence is present in the fluid any longer. Unfortunately, the GW signal from a first-order EWPT in the 
SM model plus dimension 6 operators never reaches the eLISA noise curve, not even in the most optimistic case in 
which the PT is the strongest a ~ 2.27, and the friction is the largest, 77 = 1. Note that the PT in this model cannot 
get stronger than a ~ 2.27, since higher a corresponds to smaller values of the suppression scale of the dimension 6 
operator M, for which the system becomes metastable and the Universe gets stuck in the symmetric phase [4l|. 

For small values of a, both the GW signals from bubble collision and from the MHD turbulence are sourced only 
by the bulk kinetic energy fraction k v . However, the GW spectrum from MHD turbulence dominates the signal, 
determining the global peak. The reason is simply that the GW spectrum from MHD turbulence continues to rise 
after the bubble collision peak, because it peaks at a frequency which is higher by a factor see Eq. (fl9|) . We also 
remark that as a increases, fl/H* and T* both decrease, so that the global peak of the GW spectrum is at smaller 
and smaller frequency (see Eq. (|19p and Eq. (jTTJ)) , and progressively gets shifted away from the region of maximal 
sensitivity of eLISA. This fact had already been remarked in Ref. [411 ], which for the first time correctly determined 
/3/if* and T* from the Euclidean three action, without treating them as parameters independent of a. This behaviour 
is not peculiar of the EWPT model under analysis here, but comes from the general shape of the three action Ss(T) 
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Holographic Phase Transition 




f (Hz) 

FIG. 3: GW spectra for the holographic PT [13, HH. The black line shows the eLISA noise curve. The solid lines are the GW 
spectra for T* = 100 GeV, and from top to bottom /3/H* — 6 and /3/H* = 15; the dashed lines are for T* = 10 4 GeV, and from 
top to bottom again fi/H* = 6 and fi/H 9 = 15. 

as a function of temperature: it is therefore to be generically expected for thermal first-order PTs. We finally remark 
that the double peak structure characteristic of the GW spectrum arising from the sum of the bubble collisions and 
MHD turbulence signals had already been pointed out in Ref. [29(: however, in the present analysis it is the peak due 
to the MHD turbulence which occurs at higher frequency, the one from bubble collisions being related to the duration 
of the PT 1 1 \i while the one from MHD turbulence to the size of the bubbles i?* . This is different from what is found 
in (2i|, where Eq. (ll"4l) was used, which states that the peak frequency for the GW spectrum from MHD turbulence 
is related to the turbulent eddy turnover time and not the inverse size of the bubbles 1/i?* as in Eqs. (fTS")) and ([T9|) . 

We have seen that the EWPT in the SM extended by dimension 6 operators is not strong enough to give rise 
to a GW signal which overcomes the eLISA noise curve. Moreover, to determine the GW spectrum by bubble 
collision, we have used the result of the numerical simulations [261 ]. which is in principle only valid in the envelope 
approximation. For a weakly first-order PT, the bubble wall velocity can be relatively small causing the bubbles 
to expand as deflagrations: for example, in the model under analysis, for a — 0.128 and a = 0.201 for r\ = 1, the 
bubble wall velocity is respectively Vb = 0.09 and Vb — 0.2. It is not clear whether the numerical simulations are 
adapted to describe such a situation, since the validity of the envelope approximation most probably breaks down. 
Therefore, as a second example of a PT, which instead can give rise to a detectable GW signal for eLISA, we consider 
the holographic PT corresponding to the stabilisation of the radion in the Randall-Sundrum model, first analysed 
by Randall and Servant 2007 [47| and then by Konstandin et al 2010 [48|. The appealing feature of this PT is its 
extreme supercooling: the PT is strongly first-order, so that a attains very large values leading to a high amplitude 
in the GW signal. The validity of the envelope approximation is not put int o q uestion in this case. The production 
of GW has been analysed already both in [43] and (H)]: in particular, Ref. [48| uses Eqs. ([To]) and (IT71) to estimate 
the GW signal. We show the corresponding GW spectrum in Fig. [3] In this case, the PT is so strong that the bubble 
propagation is always in the runaway regime [481 ] ; for reasonable values of the friction, we can therefore assume that 
the gradient energy of the Higgs field strongly dominates the energy balance of the PT, and set k ~ ~ 1 and 
k v ~ Pther/ Pvac — 0. Since there are no bulk motions of the fluid, no contribution from MHD turbulence is expected 
for such a strong PT, and the bubble wall velocity equals the speed of light Vt, = 1. Moreover, from Eq. (|16l) we see 
that for a > 10, the amplitude of the GW spectrum becomes virtually independent from a. The only parameters 
one is left with to estimate the GW signal are therefore /3/iJ* and T*. The parameter P/H* can take a wide range of 
values, depending on the rank of the gauge group in the dual CFT theory (c.f. table I in [11]). In Fig. [31 we show 
four examples of GW spectra for a choice of parameters for which the signal is much higher than the eLISA noise 
curve: respectively, T* = 100 GeV, = 6 and P/H* = 15, and T* = 10 4 GeV, P/H* = 6 and P/H* = 15. 

From the two examples above, it appears that once a model for the PT is chosen, a, P/H* and T* are in general 
not independent among each other; in particular, as a increases, P/H* and T* decrease. Indeed, this behaviour is 
due to the shape of the Euclidean three-action as a function of temperature: therefore, we can expect it to hold quite 
generally, and to occur in many different models of the PT. However, this relation among the parameters cannot be 
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FIG. 4: Contour plots in the (a ,/3/H*) parameter space for different values of T» and two values of the friction 77. eLISA is 
most sensitive to PTs occuring at temperatures of the order of 10 TeV, for which it can detect a wider range of values of a and 
P/H*. This means that if the PT takes place at high temperature, eLISA can detect its GW signal even if the PT is not very 
strong and occurs quite fast. The shadowed regions correspond to a signal to noise ratio of at least one. 



specified analytically or numerically in more precise terms without knowing the three-action in details, i.e. without 
restricting to a given model. Therefore, it is not possible to perform a model independent analysis which realistically 
accounts for the relation among a, and T*. As a consequence, in the following we have chosen to let them vary 

freely, and in Fig. 2] we present contour plots in the plane (a , 0/H*) for different values of T*. Note that Vb(a, 77) and 
n(a, rj) are evaluated from the analysis of Espinosa et al. 2010 23], so we do not make use of the Jouguet detonation 
hypothesis. The criterium we adopted in Fig. [4] is that we consider a GW spectrum 'detectable' (shadowed region) as 
soon as it overcomes the eLISA noise curve. However, the fact that for a given set of the parameters (a, /3/H*, T») 
the corresponding GW signal is detectable does not mean that there can exist a realistic EWPT model in which this 
set of parameters is allowed. 

The first thing to notice about the contours in Fig. |4] is that small values of /3/.H* are easier to detect: this is 
because the amplitude of the GW spectrum is quadratic (linear for MHD turbulence) in iJ*//3, c.f. Eqs. pl)]) and (|18p . 
However, this happens only for high enough temperature: if P/H„ is small and the temperature as well, then the 
signal peaks at too small frequency to hit the eLISA noise curve (we have seen that it is the case, for example, for the 
SM plus dimension 6 operators) . Correspondingly, if becomes too high, the signal peaks at too large frequency and 
the same value of /3/-H* is no longer detectable: therefore, the shadowed region does not grow monotonically with T*. 
Still, the detection region increases for increasing until T* ~ 1000 GeV, because of the shape of the eLISA noise 
curve: given a value for /3/H*, the amplitude of the spectrum can be too small for detection if the corresponding peak 
frequency is low; increasing the temperature does not affect the amplitude but shifts the peak to higher frequencies, 
causing it to enter the eLISA noise curve, which decreases as a function of frequency for small frequency. For the 
smaller value of the friction 77 = 0.2, the signal is less visible at small T„: this is due to the fact that, for the same 
fl/H*, the GW spectrum from MHD turbulence peaks at higher frequencies with respect to the bubble collision one. 
If the friction is small, the transfer of the vacuum energy to bulk kinetic energy is less efficient and correspondingly 
for the same value of a MHD turbulence is less important, rendering the GW spectrum less visible. This is the same 
reason for which, for 77 = 1, the curves show a peak at a ~ 1: for fixed /3/-H*, when a reaches the value of about one, 
the MHD turbulent peak enters the eLISA noise curve. But if a grows further, the amplitude of the MHD turbulent 
part of the spectrum is less and less important (k v — > 0) and for the same value of /3/H* the signal is no longer 
detectable. For rj = 0.2, the peak occurs at smaller a ~ 0.3. We also remark that the detectability of the signal stops 
increasing with a for sufficiently high values of a, since the amplitude of the GW spectrum reaches saturation because 
of the dependence as a/(l + a) in Eqs. (fT5|) and (fT5|) . Note that this effect is absent for ~ 1000 GeV, because 
for the values of j3/H* showed in the plot and for T„ ~ 1000 GeV, the peak of the GW spectrum corresponds to the 
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FIG. 5: Contour plots in the (T» ,/3/H*) parameter space for very strong PTs, a > 10, and reasonable values of the friction 
77 < 1. Again, we see that eLISA is most sensitive to PT at temperatures around T, ~ 10 TeV, where it can detect PTs 
occurring slowly (small /3/H,) with a signal to noise ratio higher than 10. 

minimum of the eLISA noise curve: even if an increase in a causes a very small increase in the amplitude, this can 
be enough for the GW s pec trum to enter the noise curve. 

From the analysis of [23jj it appears that very strong PTs have bubbles propagating as runaway solutions: the 
bubble wall velocity is simply Vf, = 1, and the energy balance is such that practically all of the vacuum energy is 
transformed into gradient energy of the Higgs field instead of bulk motions of the plasma, k — = 1, at least for 
reasonable values of the friction, say 77 < 1. Moreover, only bubble collision can act as a source of GW, and the 
GW spectra become independent on a for a > 10. For very strong PTs, in summary, one is left only with the two 
parameters fl/H* and T* to determine the GW signal (the holographic PT is an example of this). We analyse very 
strong first-order PTs a > 10 in Fig. [5l which shows contour plots in the parameter space (j3/H*, T*) for different 
values of the signal to noise (SN) ratio. The SN ratio has been evaluated by maximizing, for each frequency, the ratio 
among the value of the GW signal and the one of the eLISA noise curve at that frequency. The criterion adopted in 
Fig. 2] would correspond in this language to a SN ratio of at least one. We see that the SN grows for small values 
of /3/.H*, since this corresponds to high amplitudes of the GW signal. The signal is visible for values of /3/H* which 
first increase as a function of T*, and then decrease for T* > 1000 GeV: this is the same behaviour observed in Fig. 0] 
and it is due to the position of the peak of the GW spectrum, whether or not it falls in the eLISA frequency band. 
Moreover, it appears that the SN ratio is maximal for T* of the order of 10 TeV and small values of /3/i?*, which 
maximize the amplitude: in this case, the GW signal falls in the minimum of the eLISA noise curve. 



Cosmic strings are linear concentrations of energy that have a cosmological size and are produced at early-universe 
phase transitions in a variety of high-energy physics scenarios [59l - l6lT ] . They are relics of the spontaneous breaking 
of symmetries that can arise for instance in grand-unified and super-symmetric theories. A traditional example is 
provided by the Abrikosov-Nielsen-Olesen flux tubes associated to the breaking of a local U(l) symmetry in the 
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Abelian-Higgs model [62j. More recently, following an original idea in [63jj , it has been appreciated [64j-)66| that 



cosmic strings occur also as fundamental objects of string theory, known as cosmic super-strings, when the necessary 
stability requirements are satisfied Q . The formation of cosmic strings depends on the kind of symmetries that are 
broken in the early universe (they form whenever the vacuum group manifold is not simply connected), but not on 
the dynamics of the corresponding phase transitions: the phase transitions can be first-order or second-order and they 
can occur equally well during the thermal evolution of the universe or at the end of inflation, e.g. at the end of hybrid 
inflation for field theory cosmic strings and at the end of brane inflation for cosmic super-strings. Contrary to stable 
domain walls and monopoles, which would come to dominate the expansion of the universe and lead to cosmological 
problems, cosmic strings evolve so as to contribute only to a relatively small fraction of the total energy density of the 
universe at any time. Yet they involve concentrations of very high energy densities with relativistic velocities, which 
can have a variety of astrophysical and cosmological consequences, notably GW production. The detection of cosmic 
string signatures would provide a spectacular probe of the physics at the highest energy scales, well above the reach 
of any accelerator experiment. 

After being produced in the early universe, and after a possible early friction-dominated evolution (where the motion 
of the strings is damped by their interactions with the thermal plasma), a network of cosmic strings evolves towards 
an attractor solution called "scaling regime" . In this regime, the correlation length of long string segments divided by 
the Hubble radius, as well as their energy density divided by the total energy density of the universe, remain constant 
in time. In order to scale, the network of long strings must continuously loose a significant part of its energy, and it 
is widely believed that this is achieved through the production of closed loops - i.e. strings whose length is shorter 
than the Hubble radius. The typical size of the loops when they are produced, however, is still uncertain by many 
orders of magnitude (see below). Loops are produced when long string segments cross themselves and reconnect. The 
loops may then reconnect to the long string network or fragment into smaller loops, but eventually a population of 
isolated loops remains. These loops oscillate relativistically under the effect of their tension and shrink by emitting 
GW (together with other possible kind of radiation and particles, depending on the cosmic string model). The net 
result is that a significant fraction of the energy density of cosmic strings is continuously converted into GW, which 
then redshift until the present epoch. The resulting spectrum covers a wide frequency range and can be relevant for 
different kinds of GW probes, in particular pulsar timing observations, eLISA and ground-based interferometers. 

The observational signatures of cosmic strings vary in particular with two fundamental parameters that depend 
on the underlying high-energy physics scenario: their tension fi (which is equal to their energy per unit length) and 
their reconnection probability p. The gravitational effects of cosmic strings depend essentially on the dimensionless 
parameter G/x, where G is Newton constant. The reconnection probability is often p = 1 for field theory cosmic 
strings, but it can be smaller than unity for cosmic super-strings [67j |. This is expected to increase the abundance of 
strings (although the precise dependence is still uncertain, see below) and the chances to detect their signatures. There 
are already a number of constraints on the string tension, which depend more or less strongly on the parameters that 
describe the network evolution. For instance, constraints from the Cosmic Microwave Background and the Large Scale 
Structures require typically Gfi < few x 1CP 7 for p = 1, see e.g. [glj]. This bound is conservative in the sense that it 
does not depend on the properties of the loops. As we will see, current pulsar timing data lead already to significantly 
stronger upper bounds on the string tension in the case of large initial loop sizes. Note that the production of particles 
by cosmic strings leads also to complementary constraints on small values of the string tension, see e.g. [69) . 

The stochastic GW background from cosmic strings has been extensively studied in the 80's and in the beginning 
of the 90's [70l - l77j . see in particular Caldwell and Allen [76| for a detailed review of early works. More recently, it has 
been pointed out by Damour and Vilcnkin [78| (see also [79j]) that the GW signal from cosmic string loops includes 
strong infrequent bursts that can be looked for individually and should not be included in the calculation of the 
stationary and nearly Gaussian background. These GW bursts are produced by string configurations such as cusps 
(highly boosted pieces of loops where the string doubles back on itself) and kinks (shape-discontinuities that propagate 
along the strings at the speed of light). The sensitivity of the GW signal to different cosmic string parameters has 
been further studied in [8(| and the detection of individual GW bursts in [8l| . The stochastic background from loops 
was subsequently studied by Hogan et al [82| and Siemens et al [Hj|. The contribution of the long strings to the 
GW background has also been recently re-analyzed in [84], where it was found to be sub-dominant compared to the 
contribution of the loops except for very small initial loop sizes and in certain frequency range (e.g. at CMB scales, 
where the GW spectrum from loops is negligible). 

In this Section, we will study the stochastic GW background produced by cosmic string loops. Our interest in 
this paper is on stationary and nearly Gaussian GW backgrounds, so we will not consider the detection of individual 
bursts. We will however study the effect of removing the rare bursts from the stochastic background. Our calculation 
method will be mostly similar to Siemens et al [83j |. but with some differences that we will explain as we proceed. 

Three general aspects must be modeled in order to compute the GW spectrum from cosmic strings. First, we 
need a model for the expansion history of the universe. The impact of this has been studied in a more and more 
sophisticated way in the early works. For instance, Caldwell and Allen [76| took into account the variation of the 
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number of relativistic species in the thermal plasma of the early universe as it cools, as well as the effect of the energy 
density in cosmic strings and in GW on the expansion of the universe. On the other hand, Refs. [HI, [83| used more 
recent values of the cosmological parameters, considering in particular the late-time acceleration of the expansion of 
the universe. Here we will consider both a A-CDM cosmological model and the variation of the number of relativistic 
species in the early universe (see section [ill B|) . We will also discuss the dependence of the GW spectrum on the 
cosmological history. We will neglect the contribution of cosmic strings and GW to the total energy density of the 
universe. Given the improved upper bounds on the string tension, we expect this to be a good approximation. 

A second aspect that must be modeled is the number density of loops in the universe, which depends crucially 
on the typical size of the loops produced by the long string network. This issue is still under debate, with different 
simulations [85-88] and analytical works (89M9i| obtaining different results. The initial length of a loop (i.e. before it 
shrinks by emitting GW) produced by the long string network at cosmic time t can be written as 

k = at (23) 

where the length is defined as the energy of the loop in its rest frame divided by its tension (the initial length in the 
cosmological frame is modified accordingly if the loops are produced with a large initial velocity). The first simulations 
indicated that loop production occurs at relatively small scales, close to their resolution limit (see [H3] and references 
therein). This suggested that li should be given by the size of the small scale structure on the long strings, which are 
damped by gravitational emission. This led to the estimate [§2[ a ~ TG/i where r ~ 50 measures the efficiency of 
GW emission. More recently, Ref. (93| pointed out that gravitational emission may be much less efficient in damping 
the small scale structure, leading to much smaller initial l oop sizes: a oc (Gfi) x with x > 1 89]. This was further 
investigated in [9(j, which obtained different values than [89| for the exponent \. On the other hand, the latest 
simulations now indicate that is instead given by the large scale properties of the network, leading to much larger 
initial loop sizes: a ~ 0.1 [87t [88| . Finally, it is possible that loop production occurs simultaneously at two very 
different scales, a < TGfi and a ~ 0.1, as advocated in e.g. |9l| . 

Whatever the typical scale of loop production, there is in practice a distribution of different values of a around some 
peak, but this is usually approximated by a Dirac distribution in the calculations of the GW spectrum. Caldwell and 
Allen frpt ] consider different possibilities with a < 10~ 3 . Siemens et al [83[ study both a < TG/i and a ~ 0.1. Hogan 
et al 82] consider a ~ 10~ 6 — 0.1 and a case where loops are produced simultaneously at a ~ 0.1 and a ~ 10~ 5 . For 
the same values of the other parameters, the case a ~ 0.1 leads to a much larger amplitude of the GW spectrum than 
the case a < TG/j,, because the loops are long-lived (compared to the Hubble time) in the first case and short-lived in 
the second one. We will consider both cases separately below. For a < TGfj,, we will use the same loop distribution 
as in [HI (see section Ull C|) . For a ~ 0.1, our loop distribution will differ from the one used in 83] in order to account 
for the smooth transition between the radiation and matter eras and for the effects of the variation of the number of 
relativistic species and of the late-time acceleration of the universe expansion on the loop number density (see section 
UTTPl) . 

Finally, a third aspect that must be modeled in order to compute the GW background from cosmic strings is the 
GW spectrum emitted by each individual loop. As it oscillates quasi-periodically, a loop of length I emits GW at 
frequencies f em = 2n/l, where n = 1,2,3,... denotes the harmonics of the loop oscillation. The total GW power 
emitted by the loop is given by 

P = Y, p nG^ 2 =TG^ 2 (24) 

n 

where P n and V are dimensionless coefficients that depend on the loop oscillation. In practice these coefficients 
can take different values for different loops, but some average values (over the loop population) are assumed in 
order to calculate the GW spectrum. It is not clear, however, what the average values for each P n would be in a 
realistic population of loops. Different works have used different assumptions, typically all the power emitted in the 
fundamental n = 1 mode or a given asymptotic behavior P n cx n~ 1 ~ q for n — > oo (extrapolated down to n = 1). 
The value of q then depends on the presence of cusps or kinks, which are responsible for the emission at n > 1. 
Caldwell and Allen 76] consider mostly q — 1/3, which was obtained analytically for loops with cusps in p72| and in 



the numerical simulations of [9J| . Ref. [77| considers also other values of q as well as the effect of a cutoff for the power 
emitted at large n. As we already mentioned above, it was later emphasized by Damour and Vilenkin p78j that cusps 
and kinks lead to strong infrequent GW bursts that should not be included in the computation of the stationary and 
nearly Gaussian background, in order not to over-estimate its amplitude. The GW spectrum from cusps and kinks 
after removing the rare bursts was then studied by Siemens et al [83| . On the other hand, Hogan et al [82| considered 
all the power to be emitted by the first few harmonics, mainly the fundamental mode n = 1. In the following, we will 
consider different models for the GW spectrum emitted by each individual loop, and find that they only lead to small 
differences in the GW background today. We will also find that, in the case of large initial loop sizes, removing the 
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rare burst has practically no effect on the present-day GW spectrum, at least when the number of cusps and kinks 
per loop oscillation period is of order unity. 

In addition to the three general aspects discussed above, and depending on the kind of cosmic strings that one 
considers, other properties may affect the GW signal. For instance, global strings (associated to the breaking of a 
global symmetry) may radiate mainly Goldstone bosons, in addition to GW. The amplitude of the GW signal from 
global strings is then reduced, as estimated in e.g. [77J. If different kinds of strings are present and interact with each 
other (as for instance for cosmic super-string networks involving F-strings, D-strings and their bound states), the 
evolution of the network and its properties may also be more complicated. In particular, junctions between strings 
can form, leading to a drastic increase in the number of (sharp) kinks due to the fact that each kink gives rise to 
two transmitted and one reflected kink every time it reaches a junction. The impact of this phenomenon on the 
gravitational wave spectrum has been studied in [95]. It has also been found in [96j that the GW signal from cusps 
and kinks is reduced for cosmic super-strings that can move in the extra-dimensions. On the other hand, it has been 
argued in [67( that the position of the strings in the extra dimensions correspond to worldsheet moduli that should 
be stabilized, so that the classical motion of the strings should be effectively four-dimensional. Cosmic strings may 
also appear in hybrid defects, such as strings ending on domain walls or monopoles, in which case the network is 
typically unstable. The GW signal from such defects has been studied in [97l [9|| . Finally, it has been argued in [99| 
that the gravitational clustering of loops in the galactic halo may lead to specific GW signatures. In the following, 
we consider only stable cosmic string networks, which are effectively four-dimensional and whose loops decay mostly 
into GW. We also assume that only one kind of string is present or dominate the loop number density, and we neglect 
any clustering of loops. 



B. Calculation of the GW Background 



In this sub-section, we detail the procedure that we follow to compute the GW background from cosmic string 
loops. We derive an expression for the GW spectrum using two different methods. The first one is standard [6(| and 
is convenient to obtain the overall normalization of the GW spectrum. It can be applied to a variety of models for 
the GW spectrum emitted by each individual loop. The second method is similar to the one used in [78|, [83[ and is 
convenient to remove the rare bursts from the GW background. It is based on the asymptotic waveform emitted by 
cusps and kinks. As we will see, except for the overall normalization and the removing of the rare bursts, the two 
methods are equivalent. 



1. Models for the GW Spectrum Emitted by Each Individual Loop and for the Cosmological Evolution 

As discussed above, a loop of length I emits GW at frequencies f em = 2n/l at the time of emission, with the power 
(|24|) where n = 1,2,3, ... denotes the harmonics of the loop oscillation. Approximating the discrete spectrum by a 
continuous one, we have 

dP f 

— = Gti 2 l S(f em l) with / dx S{x) = T . (25) 

One common model of GW emission is provided by 

S q ( x ) = Q(x - 2) with A q = 2 q qT (26) 

which corresponds to an asymptotic behavior P n cx n~ q " 1 at large n. The numerical pre- factor A q results from the 
normalization (|2~51 and the Heaviside function (O(x) = for x < and 1 for x > 0) results from the condition n > 1 
on the harmonics. As we will see, the contribution of cusps to the asymptotic power at large n goes as P n oc n~ 4 / 3 , 
while the contribution of kinks goes as P n oc n~ 5 / 3 . Since the contribution from cusps decreases more slowly than the 
one from kinks when n increases, one expects the total power at large n to be dominated by the cusp contribution 
if the number of cusps is not strongly suppressed compared to the number of kinks. We will then use the model 
(l26t with q = 1/3 for loops with cusps, independently of the presence of kinks. We call these "cuspy loops" in the 
following. On the other hand, we use the model ([2"u| with q = 2/3 for loops that have kinks but no cusps. We call 
these "kinky loops". Note also that the model (f2l)f may be inaccurate for the low harmonics (it assumes that the 
asymptotic behavior of the power is valid down to n = 1, which need not be the case). It will therefore be interesting 
to also calculate the GW background in another model that assumes instead that all the power is emitted by the 
fundamental mode n = 1. This corresponds to 



Sfun(x) = T6(x- 2) 



(27) 
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where S(x) is the Dirac distribution. 

Consider now the energy density dp gw emitted in GW with frequencies in the interval [f em , f em + df em ] during the 
time interval [t, t + dt] by loops whose length belongs to the interval [I, I + dl]. The number density of such loops is 
denoted as n(l,t)dl. Each of them emits a GW energy dE gw = Gp, 2 lS(f em l)dt df em . This gives 

dpg W = Gp 2 I S(f em l)n(l,t) dtdldf em . (28) 

After their production, the GW energy density redshifts as 1/a 4 and their frequency as 1/a, where a is the scale 
factor. We denote by ao the value of the scale factor today and by / = f em a(i)/ao the present-day GW frequency. 
Integrating over production times and loop lengths, the present day GW spectrum is then given by 

**--<V f *^ r*lln(l,t)s( a 4f). (29) 



d f Jt, <4 Jo ' \a{t) 

Here to denotes the present time and f* the time when GW start being produced, e.g. the end of the friction-dominated 
era. Assuming that all the loops produced at time t have an initial length Zj = at before shrinking by GW emission, 
there are only loops with length < I < a t at time t. The quantity of interest is ((6]). Using 3H 2 = 8ttG p c , where Ho 
denotes the Hubble constant today, we have 

*°~-t r'"*" 5 ^) - (3o » 

In the following, it will be more convenient to work in terms of redshift z, 1 + z = Oo/a(t). We consider a A-CDM 
cosmological model where the Hubble rate at redshift z is given by 

H(z)=H H(z) with H(z) = V^a + fiji/ (1 + z) 3 + M R G{z) (1 + z) 4 . (31) 

For the cosmological parameters, we use the values given in [lOOt : H = 100 h km/ s/Mpc with h = 0.72, n A = 0.74 
for the cosmological constant abundance today and Qm = 0.26 for the present-day abundance of cold dark matter and 
baryons. In the expression for H(z), the last term inside the square root dominates during the radiation era, when 
the energy density of the thermal plasma goes as pr{z) oc g*(z) T 4 where g*{z) is the effective numbe r of relativistic 
degrees of freedom at redshift z and T is the temperature of the photon background at that time [ 1 ll j . Entropy 
conservation gives T oc gs(z)^ 1 ^ 3 a -1 where gs(z) is the effective number of entropic degrees of freedom at redshift 
z. The energy density in radiation can then be written as 

p R (z)=p R (0)g(z)(l + z) 4 with Q(z) = 9 *[ Z ] g f 3 {0) - (32) 

The function Q{z) varies as the universe cools when species become non- relativistic and release their entropy to the 
relativistic species that are still in thermal equilibrium. The evolution of g* and gs with the photon temperature 



for the particle content of the Standard Model is given in e.g. [101| : g* = g s = 106.75 when all the Standard Model 
particles are relativistic, <7» = g s = 10.75 in between the QCD phase transition and electron-positron annihilation, 
and g» = 3.36, gs — 3.91 after electron-positron annihilation and neutrinos decoupling. The neutrinos can be treated 
as if they were massless because they are relativistic during the radiation era and give a negligible contribution to the 
total energy density during the matter era. The Hubble rate during the radiation era is given by pip with an effective 
abundance of total radiation today (photons and neutrinos) given by fl R = p R (0)/p c ~ 8 x 10~ 5 for To = 2.725 A'. 
The function G{z) varies mostly around T — 200 keV (electron-positron annihilation) and T = 200 MeV (QCD phase 
transition). We model it as a piecewise constant function which is discontinuous at these two moments of time. This 
approximation gives 

( 1 for z < 10 9 

Q(z) = I 0.83 for 10 9 < z < 2 x 10 12 (33) 

[ 0.39 for z > 2 x 10 12 . 

The variation of <?* and gs will introduce features in the present-day GW spectrum and reduce its amplitude at 
high frequencies. We consider this effect in a minimal way, based solely on the particle content of the Standard 
Model. Additional entropy release in the thermal plasma at high temperature, e.g. when other particles become non- 
relativistic, would further reduce the GW spectrum at high frequencies. This, however, is highly model-dependent. 
We will come back later to the dependence of the GW spectrum from cosmic strings on the thermal history of the 
universe. 
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Integrating over redshift instead of time in Eq. (|30|) . we have 

8tt h 2 ,„ n9 . [ z * dz f at{z) 



where 



/ x / , f + °° dz' , , 

^ = -W with ^ = 1 ^ 



For the loop spectrum (|26|) . this gives 



2. Background and Rare Bursts from Cusps and Kinks 

We now consider another method to calculate the GW spectrum, which is based on the asymptotic waveform 
emitted by cusps and kinks. The logarithmic Fourier transform h(f) = \ f\ J h(t) e l2 ' K ^ t of the strain h(i) produced by 
a cusp or a kink on a loop of length I can be written as |78j 



GfJLl 
r {femlY 



h q (fem,l,r) = n« (37) 



where q — 1/3 for a cusp and q = 2/3 for a kink, / em is the GW frequency at the time of emission and r is the 
physical distance from the source. Eq. (|37|) is derived in [78[ for loops that are relatively smooth (i.e. loops that can 
be characterized by a single length scale I), up to factors of order unity that depend on the shape of the cusp or the 
kink, and in the limit n = f em l/2 S> 1 (cusps and kinks contribute to the harmonics n 3> 1 of the loop oscillation). In 
this limit, the GW signal in the time domain corresponds to a burst with a characteristic time scale At ~ l// em <^ I- 
The GW emitted from a cusp is beamed into a cone of opening solid angle d£l q w irO^ , where 

O m (fem, = 77 TyfT^ • ( 38 ) 

(Jem'') 

On the other hand, a kink emits GW in a thick fan of directions, within a solid angle dfl q f=s 2Tr8 m . Although derived 
in the limit n = f em l/2 3> 1, a common assumption is that these equations remain approximately valid down to the 
fundamental mode n — 1, up to numerical factors of order unity that depend on the particular motion of the loops. 
One can then consider an approximation where the signal is modeled by Eqs. (|37l I38[) for all the frequencies with 
n = feml/2 > 1. We continue to call the signal in this approximation a "burst", even if only its high-frequency part 
is burst-like. 

A GW emitted at redshift z travels over a physical distance 

- : i with <P r (z)= (39) 



#o ' Jo U{z' 

until today and has a present-day frequency / = f em /(l + z). Its amplitude today can then be written as 

h q (.f,z,l)= V7WTZ ^ with f l ^ + z )> 2 (4°) 
<j> r (z) [fl{l + z)) H 

where the last condition results from the condition n > 1 on the harmonics. An important quantity is the rate 
d 2 N/dzdl of GW bursts observed today at the frequency / and that were emitted between redshifts z and z + dz 
by loops of lengths between I and I + dl. There are n(l,z)dldV of such loops in this redshift interval, where dV = 
4irr 2 dr/(l + z) 3 = <i'K(j) 2 dz/(H 3 (l + z) 3 T-L(z)). Each of these loops emits 2N q /l bursts per unit emission time, i.e. 
2N q /(l(l + z)) bursts per unit observation time, where N q is the average number of cusps per loop oscillation period 
T = 1/2 (for q — 1/3) or of kinks propagating on the loops (for q = 2/3) . Only a fraction c?f2 9 /(47r) of the bursts is 
emitted towards a given observer, where dtt q was given above. This gives 

^(/, /) = ^m^hM i e {f i(i + z)-2). (4i) 

dzdl U > ' > H$lH(z)(l + z)UMl + z)) q 
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The present-day spectrum of the GW background resulting from the superposition of many bursts can then be 
calculated as [zl El H3 

h?n gw = ^^f Jdz Jdl h 2 q (f, z, I) ^(/, z, I) . (42) 

As discussed above, Ref. [78| pointed out that the strongest bursts may be rare enough to be observed individually and 
that these strong infrequent bursts should not be included in the computation of the stationary and nearly Gaussian 
background h 2 fl gw (the background resulting from the superposition of many bursts is approximately Gaussian by 
the central limit theorem) 7 . Thus only the bursts with an amplitude h q smaller than some threshold hb are effectively 
superimposed and the domain of integration in (|4"2"j) must be restricted to these bursts only. Following [78|, |83| , we 
define the threshold hb(f) at a given frequency / (e.g. the characteristic GW frequency probed by a given observation) 
by the condition that the rate of occurrence of the bursts with amplitude h q > hb is equal to /: 

dz J dl z, I) 8 (h g (f, I, z) - h b (f)) = f . (43) 

The integration in (l4"2l is then restricted to the domain where h q (f, I, z) < hb(f). 

Eq. (|42l) can be developed by using Eqs. (|40l l4Tj) . Except for the condition h q (f,l,z) < hb(f), the result then 
reduces to Eq. (|36f with A q replaced by another pre-factor that depends on the average number N q of cusps or kinks 
per loop oscillation period. These two pre-factors are equal for an effective number of cusps or kinks given by 

1 q r 
Sir* 

In that case, the GW spectrum is given by 

87T h 2 f z ' d? rat(z) n \ 

* ^ = YWo A q (G,r I H[z){l + zf+q I dl-^le (fl(l + *) - 2) (k b{ f) - h q (f, I, z)) (45) 

where the last factor removes the rare bursts with h q > hb- 

The models (f2"fj|) and (|40p thus lead to the same GW spectrum when the condition (|44p is satisfied. In the following, 
we adopt this normalization for the effective number of cusps or kinks in terms of the factor T that characterizes the 
total power emitted in GW by each individual loop, see (|24"1) . Note that the parameter T enters also in the expression 
for the loop number density n(l, z) (see below), because it determines the lifetime of each loop. For consistency, the 
same value of T should be used in the total power emitted by each loop and in the loop number density. In other words, 
if the number of cusps or kinks varies, one must take into account the resulting effect on the loop number density. 
Otherwise, one could naively conclude for instance that increasing the number of cusps or kinks would increase the 
amplitude of the GW background today, because it increases the coefficient T in the total power emitted by each 
loop. However, this conclusion is incorrect because increasing T also decreases the loop number density. As discussed 
below Eq. (|26l) . we use q = 1/3 for loops with cusps (independently of the presence of kinks) and q = 2/3 for loops 
with kinks and no cusp. In the first (resp. second) case, variations of the average number of cusps (resp. kinks) can 
then be considered by varying the parameter T. For instance, for the value T = 50 that is often considered in the 
literature, Eq. ((44)) gives N q 2 — 3. Of course Eq. (|44]) is specific to the model (f26l |40|) for the GW spectrum emitted 
by each individual loop. The number of cusps or kinks on the loops could also be suppressed without decreasing the 
parameter T, by increasing the power emitted by the low harmonics. In this case, however, the model (f2"6l |4"0")) is not 
accurate. It will therefore be interesting to also calculate the GW spectrum with the model ||2"T|) where all the power 
is emitted in the fundamental mode. 

There are few differences in Eq. (|4~5]l compared to the corresponding expression used by Siemens et al [83[ . Ref. [83[ 
considers one cusp and one kink per loop oscillation period, but uses an independent value of T in the loop number 
density. For T = 50, this decreases the GW amplitude compared to the one we will obtain below by a factor of 2 to 
3. Ref. [83[ also considers all the frequencies with fl(l + z) > 1, while we impose fl(l + z) > 2 that follows from 
the condition n > 1 on the harmonics. This slightly increases the GW amplitude obtained in~[83j compared to the 
one we will obtain below, in particular at high frequencies where the GW spectrum is dominated by the contribution 
of the low harmonics. The second paper in Ref. [HI also uses an expression analog to Eq. (|4"5]) to calculate the 



7 The transition regime between the infrequent bursts and the Gaussian background has been recently studied in 102] . 
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contributions from cusps with q = 1/3 and from kinks with q = 2/3, and then sum the two contributions. However, 
both contributions involve also the low harmonics, whose effect on the GW spectrum is then effectively counted twice 
in the model (|26l40j) . Here we consider instead that the parameter q depends on the kind of loops that would be most 
abundant in a realistic population of loops. As discussed below Eq. (j2"6")l , we then use q — 1/3 to model cu spy loops 
and q = 2/3 to model kinky ones, but we do not sum the two results. Other differences with respect to 83] come 
from our different cosmological model ([3T1 l33j) with the variation of the number of relativistic species, and from the 
different loop number density that we use in the case of large initial loop sizes (see Section UlI D|) . 

To proceed further, we need to specify the loop number density n(l,z). As discussed below Eq. (1231) . this depends 
on the typical size of the loops produced by the long string network. In the following, we consider separately the cases 
of small initial loop sizes (a < TGn) and of large initial loop sizes (a ~ 0.1). 



C. Small Initial Loop Size 

1. Loop Number Density 

As discussed in Section IIII Al the typical size of the loops when they are produced by the long string network is 
still uncertain. If their size is set by the small-scale structure on the long strings, one expects a < TG/^, where a is 
defined in Eq. ()23|) . When a loop of initial length I = aU is produced at time U , its instantaneous energy fil decreases 
by GW emission at the rate given by Eq. (pM|) . Its length at a later time t > t i is then given by 

l = ati-TGfi(t-ti). (46) 

The lifetime of the loop is obtained by setting I = in this expression, At ~ ati/(TGfi). Thus for a < TGfi, the loops 
decay in less than one Hubble time after their production. 

In that case, one expects the approximation where all the loops at time t have the same length I = at to be relatively 
accurate. We then use the same loop distribution as in [|o|, |83| 

n ^ )= mii^)i (47) 

pi Gfit ti (z) 

where p is the reconnection probability, t(z) is given in Eq. (|35|) and 

C (») = C M + {Cr - ° m) Z (48) 

Z "T Z & q 

with Cr = 10, Cm = 1 and z eq ~ CIm/^r is the redshift of radiation and matter equality. 

This expression for the loop number density is obtained by matching the energy density lost by the scaling network 
of long strings in the radiation and matter eras to the energy density produced in loops (assumed to have all the same 
size) and by taking into account the (short) lifetime of the loops, see e.g. [8(|. One then finds that the loop number 
density is about ten times larger in the radiation era than in the matter era, which is modeled by the function C(z). 
For small loops, we assume that the late-time acceleration of the universe expansion and the variation of the number 
of relativistic species in the early universe do not significantly affect the efficiency of loop production, so that Eq. (|4"7|) 
remains a good approximation 8 . 

Note that the dependence of the loop number density on the r econ nection probability in Eq. (|47l) is still uncer tain . 



Simple arguments suggest n(l, z) cx l/p 13 with /3 = 1 [80| (see also [l03j |). but the the values /3 = 2 [65( and ft = 0.6 104 1 
have also been obtained. In Eq. ([47]) and in the following, we focus on the case n(l, z) cx l/p for definiteness. However, 
the results can be directly generalized to the case where n(l,z) cx l/f(p) with an arbitrary function f(p) by simply 
replacing p by f(p) in all the expressions below. 

As discussed below Eq. (1231) , in the case of small initial loop length, it is still unclear what this initial length should 
actually be. Following [8fJ,[83|], we use the parameter 



The evolution of the long string network and the efficiency of loop production are in principle affected by the late-time acceleration and 
by the variation of , but to the best of our knowledge this has not been studied in numerical simulations. The variation of g* results 
in transient deviations from an exactly radiation-dominated expansion, which can only affect the loop number density for small loops 
in a small and transient way. 
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to parametrize this uncertainty. We will mainly treat e as a free but constant-in-time parameter, as in [8(| IH|. The 
latter assumption may however be inaccurate, as the results of [1^, [9(| indicate instead that e takes different values 
during the radiation and matter eras. We will therefore also discuss below the GW spectrum computed in such a 
case. 



2. Calculation of the GW Spectrum 

The approximation that all the loops have the same length I = at(z) at redshift z significantly simplifies the 
expressions in Section [til Bl The RHS of Eq. (l40l) decreases when z increases, so the condition h q (f, I, z) < hb(f) that 
removes the rare bursts in Eq. (|45|) corresponds to a lower bound on the redshift, z > The value of Zf>(/) is 

found from the condition (j43|) . Using (|41[) with ([47]) . this gives 



dz 



C{z)4> 2 r {z) 



pTGfi ( af 



2-q 



(50) 



where N q is given in (l44l) . The bursts produced at redshifts z < Zb(f) are too rare to contribute to h 2 tt gw (f). On the 
other hand, the condition fl(l + z) > 2 on the harmonics of the loop oscillations in Eq. (|4"5|) leads to an upper bound 
on the redshift, z < Z\{f) where 



(l + zi(/)) <t> t (zi(f)) 



2H Q 



(51) 



Compared to a given present-day frequency /, the GW produced at z > Z\{f) contribute only to the GW spectrum 
at higher frequencies. For Z\{f) < z*, Eq. (|4"S"|) then gives 



dz 



C(z) 



(52) 



This expression for the GW spectrum is valid for / < /* , where /» is a high-frequency cutoff defined by the condition 
z i{f*) — z *- Here is a high redshift above which the cosmological model (|31|33[) . or the loop number density (|47|). 
stops being applicable (see section HlIC4[) . 
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FIG. 6: Comparison of the GW spectra for cuspy loops (q = 4/3) and for kinky loops (q — 5/3), in the case of small initial 
loop size. In each case, we also show the spectrum that would be obtained without removing the rare bursts. The parameters 
are e = 1, Gfj, = 1(T 7 , p = 1 and T = 50. 

GW spectra computed with with Eqs. (l5Ul) - (|52p for a given set of parameters in the case of small initial loop sizes are 
shown in Fig. [6] The spectra are characterized by a nearly flat part at high frequency, corresponding to GW produced 
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during the radiation-dominated era, and a low-frequency peak with a very steep infra-red tail. The high-frequency 
"steps" in the GW spectrum are due to the variation of the number of relativistic species in the thermal plasma of 
the early universe (see section [Til C 4[) . In Fig. HI we compare the GW spectra obtained by using two different models 
for the GW spectrum emitted by each individual loops: Eq. (|2"B)) with q = 4/3 ("cuspy loops") and q = 5/3 ("kinky 
loops"). In both cases, we also show the GW spectra that would be obtained without removing the rare bursts. We 
see that, in the case of small initial loop sizes, the removing of the rare bursts is significant for cuspy loops, but it 
is barely noticeable for kinky loops 9 . Once the rare bursts are removed, the two models for the spectrum emitted 
by each loop lead to very similar GW spectra today. The differences occur mostly around the peak and around the 
high-frequency steps, but they always remain smaller than a factor of order 2 in the amplitude. 

G^=10 -7 , p=i 



e=0.3(Gf/)'' 2! 
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FIG. 7: GW spectrum for small cuspy loops with a different value of e in the matter and radiation eras: respectively e = eu = 
0.3 (Gfi) 1 ' 25 and e = en = 0.2 (Gfi) 0A4 (red, plain line). We also show for comparison the GW spectra for a constant value of 
e = £m (blue, dashed line at the right) and e = en (black, dashed line at the left). The other parameters are Gfi = 10~ 7 , p = 1 
and T = 50. 

As mentioned above, the results of [89|, 90] indicate that the parameter e (determining the size of the loops) may 
take different values in the matter and radiation eras, although the values obtained in these two references differ 
from each other. According to [9(| (second paper in that reference), the effective values of e that are relevant for the 
calculation of the GW spectrum (taking into account the high velocity of the loops in the cosmological frame) are 
€m ~ 0.3 (G/x) 125 in the matter era and €r « 0.2 (G^) - 44 in the radiation era. A GW spectrum calculated in this 
case is shown in Fig. [7J As we will see, when e decreases, the GW spectrum is shifted towards higher frequencies, 
while its overall amplitude remains unchanged. For a given value of e during the matter era, a larger value of e during 
the radiation era then shifts the almost fiat part of the spectrum towards lower frequencies, effectively cutting out an 
intermediate part of the spectrum, while the GW spectrum around the peak is practically not modified, see Fig. [7] 
The differences between the GW spectrum with e varying from en to cm and the GW spectrum with a constant 
value of e = cm occur essentially because of the high-frequency steps in the GW spectrum due to the variation of 
the number of relativistic species. However, the two spectra never differ by more than a factor of order two or so 
in the amplitude (this remains true for other values of G/i as well). If e takes different values in the radiation and 
matter eras, one can then apply to a good approximation the results obtained for a constant value of e = ej^r . We will 
therefore focus on the case where e is constant in the following. 



3. Comparison With Observations 

The GW background from cosmic strings may be accessible to different observations. We now study the regions 
of the parameter space that can be probed eLISA and other experiments that will operate by the time eLISA flies, 



This was also noticed analytically in the third paper of Ref . |95l . 
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see Fig. [8] For ground-based interferometers, we consider the current limit set by the LIGO S5 run and the expected 
sensitivity of Advanced LIGO, both taken from Q. For pulsar timing experiments, we consider the current limit 
obtained by Jenet et al 0] and the expected sensitivity of the complete Parks Pulsar Timing Array (PPTA), both 
taken from 0, Il05| . Regarding in particular the sensitivity of PPTA, we note that the back ground from cosmic strings 
may be hidden by the stochastic background from cosmic massive black hole binaries [l05| - or vice versa, but in any 
case it will probably be very difficult to identify the source using only pulsars in the case of a detection. For eLISA, 
the confusion background from galactic binaries is already removed from the sensitivity curve. 




FIG. 8: GW spectra for cuspy loops compared to observational sensitivities for different values of the cosmic string parameters, 
in the case of small initial loop sizes with T = 50. 



Since we saw above that the GW spectrum does not depend strongly on the particular model for the power emitted 
by each individual loop once the rare bursts are removed, we focus for defmiteness on the model (|26|) with q = 4/3 
("cuspy loops") in this sub-section. Fig. [8] illustrates how the GW spectra vary with cosmic string parameters and 
how they compare to observational sensitivities. The dependence of the GW spectrum on the parameters will be 
further discussed in Section fill C 41 When the string tension n decreases, the GW spectrum moves towards higher 
frequencies as / oc l/(G/i) (up to the small effect of removing the rare bursts), and its amplitude at high frequencies 
decreases as Cl gw oc G\i. Decreasing the reconnection probability p simply increases the GW amplitude as fl gw oc 1/p 
(again neglecting the removal of the rare bursts), as for the loop number density, see the discussion below Eq. (|48l) . 
When the size e of the loops decreases, the GW spectrum moves towards higher frequencies as / oc 1/e, and its overall 
amplitude remains unchanged. 

We show in Fig. [§] the regions of the parameter space that can be probed by each experiment in the case of small 
initial loop sizes. We set T = 50 and show the observational reaches in the (G/^, e)-plane for p = 1, p = 10 _1 , p = 10~ 2 
and p = 10 -3 . We use the simple criterion that the signal h 2 il gw (f) from cosmic strings is observable by a given 
experiment as soon as there is at least one frequency where the signal exceeds the sensitivity curve of that experiment. 
This criterion is accurate enough to determine the different possibilities, although the precise actual reach of different 
experiments will ultimately depend on both technological and data-analysis developments. 

The shapes of the detectable regions for each experiment can be easily understood from Fig. [8l For p = 1, which is 
expected for field theory cosmic strings, the high-frequency "flat" part of the spectrum is observable only for relatively 
large values of the string tension. For instance, the flat part of the spectrum falls always below the sensitivity of 
eLISA when G/i < I0 -8 . In that case, and provided that G/j, is not too small, only the peak of the GW spectrum can 
be observable. This requires that the peak falls in an appropriate frequency range, so it is only possible for a limited 
range of values of e . For larger tensions, the flat part is above the experimental sensitivity so the signal is detectable 
regardless of the value of e (unless of course the whole spectrum gets shifted out of the frequency band, which happens 
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FIG. 9: Regions of the parameter space that can be probed by current and future observations, in the case of cuspy loops with 
small initial size and T — 50. We slice the parameter space in the (G/x, e)-plane for p — 1 (top left), p = 1CP 1 (top right), 
p — 10~ 2 (bottom left) and p — 10~ 3 (bottom right). The darkest regions are already ruled out by current observations (either 
pulsars or LIGO S5 run) for these values of the parameters. The results are shown in the case where the loop number density 
is inversely proportional to the reconnection probability, n(l,z) oc 1/p, see the discussion above Eq. (|49[) . In the case where 
n(l,z) oc l/f(p) with an arbitrary function f(p), the same results are obtained for f(p) = 1, 10" 1 , 10" 2 and 1CT 3 . 



for very small values of e). The regions of the parameter space that can be probed by Advanced LIGO and pulsar 
experiments are determined by similar considerations, except that pulsars probe such a low frequency band that the 
flat part of the GW spectrum always starts at higher frequencies, for any values of the parameters. As p decreases, 
the amplitude of the whole signal increases and the main effect is to push the regions of detectability to the left: lower 
values of G/i now become observable. As a side effect, the observable regions are also shifted upwards: when the 
string tension decreases, the GW spectrum moves towards higher frequencies, which shifts the values of e for which 
the peak of the spectrum falls in the frequency range of a given experiment. 

We see that PPTA, eLISA and Advanced LIGO are complementary in the search for the stochastic background 
from cosmic strings. The signal can be simultaneously detected by the three experiments for sufficiently large values of 
the Gfi and e. This would enable to characterize the signal with good precision and increase our ability to distinguish 
it from other stochastic backgrounds. For smaller values of e, but still large enough string tensions, the GW spectrum 
moves away from the frequency band of pulsar observations, but it remains observable by both eLISA and the advanced 
versions of the ground-based interferometers. For very small values of e, the spectrum also moves away from the eLISA 
frequency band, so that the only hope to detect the signal relies on ground-based interferometers. Finally, for smaller 
values of the string tension, each experiment probes a different range of loop sizes: pulsars are sensitive to values of 
e that are few orders of magnitude smaller than e ~ 1, eLISA to an intermediate range of smaller values, and the 
advanced versions of the ground-based detectors to very small values e. Therefore, eLISA will also probe a region of 
the parameter space that will not have been accessible to any other GW experiment. 
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4- Analytical Calculation and Dependence on the Cosmological History 

In order to further study the dependence of the GW signal on the parameters and on the thermal history of the 
universe, it is useful to compute it analytically. This is easily done for present-day frequencies / 3> f m , where 

50 10- 8 10- 7 Hz 

/m ~T-G^^— ■ (53) 

This corresponds to the high-frequency part of the GW spectrum, which is produced during the radiation era. In 
this limit, the integral in the RHS of Eq. (|52|) is dominated by its upper limit z = Z\{f) ^> z eq , and one can use 
the asymptotic behavior of the cosmological functions H(z), 4>t{z) and <p r (z) for z z eq . For a given present-day 
frequency / 3> f m , the redshift Z\(f) given in Eq. (|5ip corresponds to the temperature 

T tf) r Gfj, ef 

GeV~5010^Ife f ° r f>>fm (M) 

in the radiation era. Most of the GW observed today at a given frequency / 3> f m were produced when the early 
universe had a temperature of the order of T(f). For instance, for T = 50 and G/i ~ 10 -8 , the GW that fall in the 
eLISA frequency range today were produced at temperatures of the order of T(f) ~ e MeV or so. Note that this 
differs from e.g. GW from first-order phase transitions, which must be produced at significantly higher temperatures 
in order to fall in the eLISA frequency range. This is because the GW from cosmic strings have a much smaller 
wavelength at the time of their production, especially in the case of small initial loop sizes. 
Performing the integral in Eq. (|52[) then gives 



h 2 n gw ~ 1.1 x 1Q- 3 C R -f ( - — ) for /»/ m (55) 



Gjj, f 100 \ 1/a 
P \9*(T(f))J 

where g*(T(f)) is the number of relativistic degrees of freedom at the temperature (|54l) for a given frequency /. This 
part of the spectrum is almost flat, except for the dependence in g*(T(f)). Note also that it neither depends on the 
parameter V nor on the parameter q - describing respectively the total power ([24)) and the spectrum ([26| emitted 
by each individual loop. Our numerical results agree very well with the analytical estimate (I55|) . For instance, the 
" steps" in the high-frequency part of the GW spectrum in Fig. [6] are due to the variation of the number of relativistic 
species as the early universe cools down. The most significant step in Fig. [6] corresponds to the QCD phase transition 
at T ~ 200 MeV. It occurs around / ~ 0.01 Hz, in agreement with Eq. ([54|) . For the particle content of the Standard 
Model, the variation of reduces the GW amplitude at high frequencies by a factor of order 2.6, so the effect is 
relatively small. Depending on the physics beyond the Standard Model, additional entropy release in the thermal 
plasma at high temperature would further reduce the GW spectrum at high frequencies. 

As mentioned above, the GW spectrum is given by (1521 only for frequencies / < /*, where Zi(/*) = and z r is a 
high redshift above which the cosmological model (|3~Tll3"3l or the loop number density (l4"T|) stops being applicable. For 
instance, z* can be the redshift at the end of the friction-dominated era or when the scaling regime for the evolution 
of the cosmic string network starts, or the redshift when the standard thermal evolution of the universe starts (e.g. 
the end of reheating after inflation), whichever is smaller. Thus in general z* and thus /* are highly model-dependent. 
However, a rather general constraint on the GW background, at least for field-theory cosmic strings, is that it should 
not involve GW produced at temperatures T > Tg, where Td ~ G 1 ^ 2 ^ is the temperature at the end of the friction- 
dominated evolution of cosmic string networks [6fJ 10 . Indeed, even if cosmic strings are produced earlier in the 
radiation era, their motion is usually highly damped by their interactions with the thermal plasma at T > T^, so 
much less GW (if any) are emitted during this epoch. From Eq. (l54l , this implies that the present-day GW spectrum 
would be cut off at the frequency /* ~ 10 n Hz/e in the case e < 1. Even for e = 1, this is a very high frequency, 
well above the frequency range of e.g. ground-based interferometers. Thus, in the case of small initial loop sizes, the 
friction-dominated epoch is irrelevant for the direct search of GW from cosmic strings. 

Let us now consider the following question: what is the maximal temperature in the early universe that could be 
probed by observations of GW backgrounds from cosmic strings? From Eq. ([54)) , we see that the maximal temperature 
is obtained for e = 1, and G/i and / as high as possible. Taking G/i ~ 10 -7 as an upper bound on the string tension 
and / ~ 10 2 Hz for the frequency range of ground-based interferometers, this gives T(f) ~ 1 TeV. Thus in the 



This expression for applies only to field-theory cosmic strings. For cosmic super-strings, the friction-dominated epoch ends at higher 
temperatures, see [kJ. In that case, the high-frequency cutoff in the GW spectrum occurs at even higher frequencies. 



33 



case of small loops, the observation of GW from cosmic strings can provide informations about the early universe at 
temperatures as high as T ~ TeV. Consider the case with e = 1 and p = 1 in more detail. In that case, the string 
tension cannot be much smaller than G/i ~ 10~ 7 if the resulting GW are to be observable by eLISA or Advanced 
LIGO. The GW that fall in the frequency range of Advanced LIGO were then produced around T(f) ~ TeV, while 
those falling in the frequency range of eLISA were produced around T(f) ~ 10 — 100 MeV, i.e. a bit before the start 
of Big Bang Nucleosynthesis (BBN). The amplitude of the GW spectrum in the frequency range of ground-based 
interferometers may then be suppressed compared to its amplitude in the eLISA frequency range if the evolution 
of the universe is non-standard at temperatures between 10 — 100 MeV and the TeV. This occurs for instance if a 
scalar held (e.g. a modulus) dominates the universe expansion, or more generally if any entropy is released in the 
thermal plasma, during that epoch. This would dilute any GW produced before that moment of time and suppress 
the amplitude of the GW spectrum in the frequency range of Advanced LIGO. This can in fact be a typic al sit uation 
in super-gravity or string theory with moduli condensates decaying just before the onset of BBN, see e.g [l06j | . This 
example illustrates the kind of information we could get about the early universe if GW from cosmic strings were 
observed, and it further stresses the complementarity between eLISA and ground-based interferometers. 



D. Large Initial Loop Size 

1. Loop Number Density 

As discussed in Section MI A[ the most recent simulations [13, HH indicate that the initial size of the loops when 
they are produced is related to the large scale properties of the long string network, with a ~ 0.1 in Eq. (]23f) . In that 
case, the loops are long-lived compared to the Hubble time and there is at any time t a distribution of loops with a 
wide range of lengths, from I = to I ~ 0.1 1. From Eq. (|46|). the loops that are present at time t have been produced 
at times U as small as ti ~ TG/it/a, which is much smaller than t for a ~ 0.1. The number density of loops of size I 
at time t in that case can be calculated as in [60] . It can be written as 



with 



I + TGut 
U = — € 

a 



a 



(57) 



where the last expression follows from (pl6")) with TG/j, <C a. The first factor in the RHS of Eq. (|56"|) comes from 
matching the energy density lost by the scaling network of long strings to the energy density produced in loops of 
length aU, with a different numerical factor C(ti) in the radiation and matter eras, as in the case of small initial 
loop sizes. The second factor in the RHS of Eq. (j5"6"]) . on the other hand, comes from the dilution of the loop number 
number density with the expansion of the universe between the time ij of their production and the time t under 
consideration. 

In Eq. (1561) and in the following, we again assume for definiteness that the loop number density varies with the 
reconnection probability p as n(l,z) oc 1/p, see the discussion above Eq. (ji§|). Again, the results can be directly 
generalized to the case where n(l, z) oc 1/ f(p) with an arbitrary function f(p) by simply replacing p by f(p) in all the 
expressions below. 

Using a(t) oc t 1 ^ 2 in the radiation era and a(t) oc t 2 ! 3 in the matter era, and assuming a sharp transition between 
the two behaviors at time t eq , Eq. (1561) reduces to the standard result given in 60] for the loop number density in 
the radiation and matter eras. Here we use the more general expression (|56[) for the loop number density, because it 
applies to a more realistic cosmological evolution. This allows us to account for the smooth transition between the 
radiation and matter eras and for the effects of the variation of the number of relativistic species and of the late-time 
acceleration of the universe expansion on the loop number density: these effects are encoded in the factor a 3 (ti) j a 3 (t) 
in Eq. (|5d| 11 . This leads to some differences in the GW spectra that we obtain compared to those of Siemens et 



The variation of g* and the late-time acceleration affect the cosmological evolution and therefore the dilution of the density of loops after 
their production: this is included in the second factor in the RHS of Eq. J 56 D - They may also affect the efficiency of loop production, 
as we already mentioned in the case of small loops, and therefore the first factor in Eq. J56}. We expect however these effects on the 
production of loops to be negligible for the GW background and we do not include them. Regarding the variation of g*, this can 
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al [83j, which used the loop number density given in [6CJ: mainly a more pronounced peak, similarly to the GW 
spectra obtained by Hogan et al [82] with another method. 

Numerically, it does not cost much to work with the variables I and z instead of i, and t. We define a function 
Zi(l,z) = Z (ti(l,z)), where the function ti(l,z) is given by (|57[) with t(z) given in p5[) . and the function Z(t) gives 
the redshift associated to any time t (i.e. Z(t) = t^ 1 (z)). The number density of loops of length / at redshift z can 
then be written as 



n(l, z) 



C(l,z)a 2 



(1 + zf 



p {l + TGfit(z)) 4 (1 + Z,(l,z)f 



(58) 



where 



C(l,z)=C M + 



(C R -C M )Zi(l,z) 



Zi(l,z) 



(59) 



with Cr = 6 and Cm — 0.48. We take these numerical values of Cr and Cm from [6C 
and Cm ~ 0.12 Cm with Cm ~ 4 (more recent simulations [H| give similar results). 



C R « 0.4 Cr with Cr « 15 



2. Calculation of the GW Spectrum 

The GW spectrum can be calculated with Eqs. (|45|) where n(Z, z) is given by Eq. (|58|) . However, in order to impose 
the condition /i g (/, /, z) < hb(f) that removes the rare bursts, it is more convenient to work with the variables h q and 
z instead of I and z. We do that as Siemens et al 83], using Eq. ([4TJ]) to express I in terms of h q and z, and changing 
the variables (l,z) into (h q ,z) in the double integrals (|4"2l |4"3"]) . The value of hb(f) is then found from the condition 
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and used as an upper limit in the calculation of the GW spectrum 
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only lead to small and transient modifications of loop production, as in the case of small initial loop sizes. Regarding the late-time 
acceleration of the universe expansion, its effect on loop production can only affect the number density of large loops with Z ~ 0-li 
today. These loops are negligible for the calculation of the GW spectrum. 
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FIG. 10: Present-day GW spectra obtained with three different modefs for the GW spectrum emitted by each individual loop, 
in the case of large initial loop size with a = 0.1, G/i = 10 -11 , p = 1 and T = 50. We compare the loop emission models (|26p 
with g = 4/3 (cuspy loops), ((26} with q — 5/3 (kinky loops) and (|27[) (loops emitting only in the fundamental mode). The 
effect of removing the rare bursts is so small that it cannot be seen on the plot. 

GW spectra for a given set of parameters in the case of large initial loop sizes are shown in Fig. QJJ] The spectra 
are again characterized by a nearly flat part at high frequencies and a peak at lower frequencies. However, compared 
to the case of small initial loop sizes, the peak of the GW spectrum is less pronounced and its infrared tail is less 
steep. The latter goes as h 2 £l gw oc f 3 ' 2 and is emitted at redshifts z ~ 1 by loops that were produced during the 
radiation era 12 . The high-frequency "steps" are again due to the variation of in the early universe. In Fig. 1101 
we compare the GW spectra for large loops obtained for three different models of the GW spectrum emitted by each 
individual loop: Eq. (gSJ) with q = 4/3 ("cuspy loops"), Eq. (gBJ with q = 5/3 ("kinky loops") and Eq. (gT} (all the 
power emitted in the fundamental mode). In the first two cases, we also compared the GW spectra computed with 
and without removing the rare burst. We found no difference at all in the frequency range shown in Fig. 1101 The 
effect of removing the rare bursts starts to be noticeable (but remains small) only at lower frequencies, where the 
amplitude of the GW spectrum has already dropped significantly. Therefore, removing the rare bursts has practically 
no effect in the case of large loops. We also see in Fig. QJJ that the three models of GW emission by each loop lead 
to very similar GW spectra today. As in the case of small initial loop sizes, the differences occur mostly around the 
peak and the steps, but they remain relatively small, and the amplitude at very-high frequencies is exactly the same 
in the three cases (see Section IIIID4[) . 

3. Comparison with Observations 

We now discuss the regions of the parameters that can be probed by observations in the case of large initial loop 
size. We consider the same experiments as in Section IIII C 31 and we focus again on cuspy loops since we saw above 
that other models for the GW spectrum emitted by each loop lead to similar results. Fig. [TT1 illustrates how the GW 
spectrum varies with the parameters (see Section [ill D 41 for more details) and how it compares with the observational 
sensitivities. When G/i decreases, the spectrum moves towards higher frequencies as / oc l/(G/i), as in the case 
of small loops, but its amplitude now decreases only as Cl gw oc y/G/j, at high frequencies. As a consequence, each 
experiment can probe smaller values of the string tension in the case of large initial loop sizes. For a reconnection 
probability p = 1, eLISA will be able to reach string tensions as small as G/i ~ 10~ 13 . Inside this range of accessible 
tensions, the GW spectrum is always nearly flat in the eLISA frequency band. On the other hand, for p = 10~ 3 , eLISA 
can reach much smaller values of the string tension, down to G/i ~ 10~ 17 . In this case, eLISA probes the power-law 
infrared tail of the GW background, so the slope of the spectrum is very different, see Fig. [TT] The infrared tail of 
the GW spectrum goes as H. gw oc / 3 / 2 (see Section fill D 2jl and it can only be observed by eLISA for reconnection 
probabilities smaller than unity, more precisely for p < 10~ 2 or so. The observation of this part of the spectrum could 
therefore provide a signature of cosmic super-strings, as opposed to field-theory cosmic strings with p = 1. 

Fig. II 21 shows the regions of the parameter space in the (G/i, p)-planc that can be probed by each experiment - since 
the GW amplitude increases when G/i increases and when p decreases, the detectable regions are always located to 
the right of the curve corresponding to each experiment. In order to understand the shape of the observable regions, it 
is useful to decompose the GW spectrum into a nearly flat high-frequency part and a power-law infrared tail, ignoring 
the weak peak. The border of the region probed by eLISA exhibits roughly two different regimes. For p ~ 0.1 — 1 or 
so, the GW spectrum is nearly flat in the frequency band of interest for the the smallest value of the string tension 
that can be observed. On the other hand, for smaller values of p, eLISA probes the power-law infrared tail of the 
GW spectrum for the smallest value of G/i that it can observe, see e.g. the spectrum for G/i = 10~ 17 and p = 10~ 3 
in Fig. 1111 Since that part of the spectrum varies more strongly with G/i - it varies as Q gw oc (G/i) 2 - the border of 
the observable region in Fig. [T^] becomes steeper for small values of p. A similar argument explains the shape of the 
region accessible to current pulsar observations. On the other hand, in all the range of values of p shown in Fig. 1121 
Advanced LIGO and PPTA always probe, respectively the nearly flat part of the GW spectrum and its power-law 
infrared tail for the smallest value of the string tension they can detect, so that the border of the observable region 
for each of these experiments is essentially a straight line. 



12 At even smaller frequencies, not shown in Fig. 1101 the spectrum goes as h?Q, gw oc /. That part of the spectrum is produced at z ~ 1 
by loops that were produced during the matter era. 
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FIG. 12: Regions of the parameter space in the (Gfj,, p)-plane that can be probed by current and future observations, in the 
case of cuspy loops and large initial size with a = 0.1 and F = 50. The results are shown in the case where the loop number 
density is inversely proportional to the reconnection probability, n(l,z) oc 1/p, see the discussion above Eq. (|49p . In the case 
where n(l, z) oc l//(p) with an arbitrary function f(p), p should be replaced by f(p) on the y-axis of the plot. 



Interestingly, the GW signal can be simultaneously observable by PPTA, eLISA and Advanced LIGO for a wide 
range of string tensions, down to Gfi ~ 10~ 12 . For smaller values of G/i, the ability to observe the signal with PPTA 
quickly decreases, but there remains a wide region of the parameter space that can be probed by both eLISA and 
Advanced LIGO. Again, a joint detection would improve the characterization of the signal and strengthen the case for 
a background produced by cosmic strings. For very small values of p, Advanced LIGO can reach the smallest values of 
the string tension, because for these values the GW spectrum is significantly shifted towards higher frequencies. On 
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the other hand, eLISA will be able to explore a new region of the parameter space for p ~ 1CP 2 — 1, gaining roughly 
one order of magnitude compared to the smallest value of the string tension that can be probed by Advanced LIGO 
for these values of the reconnection probability. 



4- Analytical Calculation and Dependence on the Cosmological History 



As in the case of small initial loop sizes, the high-frequency part of the GW spectrum was produced during the 
radiation era and can be calculated analytically, see e.g. [601 ] . For a given frequency today / 3> f m , where 

/™~^^10- 9 Hz, (62) 

the double integral in ([BTj) is dominated by the vicinity of the point (z, I) defined by the conditions (1 + z)fl = 2 and 
I = TGfJLt(z). The solution z(f) of these conditions satisfies z(f) ^> z eq for / 3> f m . Thus most of the GW observed 
today at a frequency / 3> f m were produced by the low harmonics n = fl(l + z)/2 ~ 1 of the loops with length 
around I ~ TG/i t(z), which are the most abundant ones and are going to die by GW emission in less than one Hubble 
time. These GW were produced at redshifts z such that t(z) z ~ l/(fTGfi), which corresponds to temperatures of 
the order of 

GeV~50lF*Ife f ° r f>>fm (63) 

in the radiation era. This is the same as Eq. for small loops, but with e = 1. Note however that the estimate (|6"3")l 
for large loops is less accurate than the estimate (|54p for small loops, in the sense that in the first case the spectrum 
observed today at frequency / is dominated by GW that were emitted in a relatively wider range of times around the 
time t(f) when the temperature was T(f). Nevertheless, Eq. (j6"3"|) still provides a useful order-of- magnitude estimate. 
The most abundant loops (those of length I ~ TGfit) at time t(f) when the temperature is T(f) were produced at 
times U(f) ~ TGfit(f)/a, when the temperature was 

To summarize: most of the GW observed today at a given frequency / 3> f m were emitted in the early universe 
when the temperature was of the order of (|6"5|) . by loops that were produced from the long string network when the 
temperature was of the order of (|64p. For large loops, the GW spectrum observed at frequency / is sensitive to the 
thermal history of the universe up to temperatures Tj(/), instead of only T(f) in the case of small loops. 

The GW spectrum at / 3> f m can be calculated analytically as in e.g. [60, Hj]. Taking into account the variation 
of the number of relativistic degrees of freedom, we obtain 

9*{T(f)) V /4 , , 
9*( T i(f))J 



h 2 n gw ~ 7.3 x 10" 

where g*(T(f)) and £f*(Tj(/)) are the numbers of relativistic degrees of freedom at the temperatures (|63)l and (f64|) . 
respectively, and we used the values of the cosmological parameters given below Eq. (1311) . This result applies for 
any value of a > TGfi (large initial loop size). Again, the spectrum is flat except for the frequency-dependence in 
g*(T(f)) and g*{Ti(f)). The last factor in the RHS of Eq. comes from the effect of the variation of g* on the 
loop number density. It reduces to one for g*(T(f)) — g*(Ti(f)), i.e. when g* does not vary in between the moment 
where GW are emitted by a given population of loops and the moment where these loops were born. It is interesting 
to note that Eq. (|65p is independent of the particular model that is chosen to describe the GW spectrum emitted by 
each individual loop: it applies to the model ([2^1) for any value of q > 1, as well as to the model (f2"T)l where all the 
power is emitted by the fundamental harmonic. However, contrary to the result (|55[) for small loops, Eq. (1651) for large 
loops now depends on the parameter V that describes the total power emitted by each loop. Specifically, increasing 
r decreases the amplitude of the GW spectrum because, although it increases the power emitted by each individual 
loop, it also decreases the lifetime of the loops and therefore their number density. Note also that Vt gw oc \JG[i for 
large loops, while Q gw oc Gfj, for small loops. 

Again, the analytical estimate (|65|) agrees very accurately with our numerical results for the high-frequency part of 
the GW spectrum, except around the location of the steps due to the variation of the number of relativistic species. 
Compared to the case of small initial loop sizes, these steps are smoother and extend over a wider frequency range, 
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see e.g. Fig. [lOj In particular, a feature in the cosmological evolution at a given temperature T now affects the GW 
spectrum in an interval around two different frequencies, those given by Eqs. (j6"3")l and As in the case of small 

loops, the variation of leads to a total reduction of the GW amplitude at high frequencies by a factor of order 2.6 
for the particle content of the Standard Model. However, when this reduction occurs at frequencies smaller than the 
ones where a given experiment operates, it now increases by almost an order of magnitude the minimal value of the 
string tension that this experiment can probe, because £l gw oc at high-frequency for large loops. This applies 

mostly to Advanced LIGO, which operates at higher frequencies. 

As in the case of small loops, we can determine the maximal temperature in the early universe that can be probed 
by observations of GW today at a given frequency /. For large loops, this is not given by the temperature T(f) in 
Eq. (|63"1) when most of these GW were emitted, but by the temperature T{(f) in Eq. (|64|) when most of the loops that 
emit these GW were produced from the long string network. For instance, if there is a significant release of entropy in 
the thermal plasma between the temperatures T{f) and Ti(f), the number density of the most abundant loops when 
the temperature is T(/) would be suppressed, which would in turn reduce the amplitude of the GW emitted at that 
moment of time. Taking / ~ I0 2 Hz for ground-based experiments and G\i < 10 -8 as required by current pulsars 
observations, Eqs.(|63l 164)) with a ~ 0.1 and r ~ 50 then give Tj(/) ~ few x 10 4 GeV for the maximal temperature 
that can be probed by observations of the GW background from cosmic strings. Again, if the standard thermal history 
of the universe starts only at a lower temperature, the GW spectrum can be suppressed in the frequency range of 
Advanced LIGO. 

Let us also check that the friction-dominated epoch for field-theory cosmic strings does not affect the regions of 
the parameter space that can be probed by each experiment in Fig. Q21 This requires that Tj(/) < Td ~ G 1 / 2 /!. 
The strongest constraint is obtained for ground-based interferometers: G/i > 10 -21 for / ~ I0 2 Hz, a ~ 0.1 and 
r ~ 50. This is actually quite close to the minimal value of the string tension G/i ~ 10 -19 shown in Fig. Q21 but 
this value can only be reached by Advanced LIGO for p ~ 10 -3 . Such a small reconnection probability is more 
typical of cosmic super-strings, for which the friction-dominated epoch can end at much higher temperatures [69| . 
For field-theory cosmic strings with p = 1, Advanced LIGO can only reach G/i ~ 10~ 12 or so. We therefore conclude 
that the friction-dominated epoch does not affect the regions of the parameter space that can be probed by each 
experiment in Fig. 1121 




FIG. 13: Examples of GW spectra from large and small loops, together with the redshift and the temperature when most of 
the GW contributing to different parts of the spectrum were produced. Most of the loops emitting these GW were themselves 
produced much earlier from the string netowk in the case of large initial loop sizes, see the main text for details. 



As a final illustration, we show in Fig. Q2] examples of GW spectra in the cases of large and small loops, together 
with the redshift and the temperature when most of the GW contributing to different parts of the spectrum were 
produced. At high frequencies, this corresponds to the temperatures (|54|) for small loops and (|63[) for large loops. In 
the case of small loops, the GW produced at redshifts z ~ 1 contribute mostly to the region of the spectrum around 
the peak. In the case of large loops, the low-frequency tail of the GW spectrum is produced at z ~ 1, as discussed in 
section IHTTJ21 
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IV. OTHER COSMOLOGICAL SOURCES 

In addition to first-order phase transitions and cosmic strings, many other cosmological sources of GW have been 
discussed in the literature. We review several of them in this Section, with a particular attention to the predictions 
in the eLISA frequency band. Other reviews with more details about some of these sources are given in 2-5]. We 
consider GW sources that operate during inflation (section ftV Ap . just after inflation (section liVB[) and during the 
subsequent thermal evolution (section lIV Cp . Of course, since we have currently no direct probe of the universe at 
energies relevant for the eLISA frequency band, one must also keep in mind the possibility of GW produced by as yet 
unforeseen cosmological sources. 

A. Gravitational Waves from Inflation and Extensions Thereof 

One GW background that is expected with high confidence is the one generated during inflation. What is much 
more uncertain is its amplitude, which varies as h 2 il gw oc G 2 Ef n& with the (unknown) ener gy scale -Ei n fl of inflation. 
During inflation, quantum fluctuations of the graviton field are parametrica lly am plified [l07j into tensor perturbations 
at super- Hubble scales by the quasi-exponential expansion of the universe [1081 ) . The amplitude of each Fourier mode 
is directly proportional to the Hubble rate at horizon exit and remains constant during all the time spent by the mode 
outside the Hubble radius. When the tensor modes re-enter the Hubble radius during the post-inflationary evolution, 
they become standard gravitational waves whose energy density is diluted as radiation with the expansion of the 
universe. The resulting spectrum today covers a very wide range of scales, from the Hubble scale just after inflation 
(corresponding to modes produced at the end of inflation) to the present-day one (corresponding to modes produced 
earlier during inflation and that only re-entered the Hubble radius today) . The amplitude of the GW spectrum today 
at a given frequency is directly proportional to the almost constant (but slowly decreasing) energy density when the 
corresponding mode left the Hubble radius during inflation. All the modes that re-entered the Hubble radius during 
the radiation-dominated era are subsequently diluted by the same amount compared to the background energy density 
at that time and lead to an almost flat (but slightly red) spectrum. After radiation and matter equality, the later a 
mode re-enters the Hubble radius, the less it is diluted compared to the background energy density. As a consequence, 
the present-day spectrum increases when the frequency decreases as h 2 fl gw oc f~ 2 (again up to a small correction due 
to the slow decrease of the energy density during inflation) in the frequency range 10 _18 Hz < / < 10 _16 Hz. 

If inflation occurs at sufficiently high energy scales (i.e. not much below the GUT scale, E ln f[ ~ 10 16 GeV), the 
resulting GW in the frequency ran ge / ~ 10~ 18 — 10~ 16 Hz may be indirectly detected in the future through their 
effect on the CMB polarization [l09j j . Except in some circumstances to be discussed below, the current CMB bound on 
inflationary GW leads already to the upper bound h 2 VL gvl <^ 1CT 15 on their amplitude at freq uencies accessible by direct 
detection experiments (the precise bound depends on the inflationary model, see e.g. [11 01 ) ^ . This is much belo w the 
reach of eLISA and ground-based experiments. Post-eLISA missions in space, like the Big Bang Observer (BBO) flTT| 
and the Deci- Hertz Interferometer Gravitational- wave Observer (DECIGO) |l!2| |. might reach the required sensitivity 
in a frequency band around / ~ 0.1 — 1Hz, where the astrophysical foreground is expected to be easier to remove. 

A couple of ways to increase the inflationary GW spectrum in the direct detection frequency range, while satisfying 
the CMB constraint at low frequencies, have been discussed in the literature. One possibility is that the early universe 
was dominated during some epoch after inflation by a fluid stiffer than radiation, i.e. with an equatio n of state w > 1/3, 
e.g. a scalar field with a potential energy that is negligible compared to its kinetic energy (w w 1) jll3l Ill4 |. During 
such an epoch, the background energy density is diluted by the expansion of the universe as a~ v with v > 4, while the 
GW energy density still redshifts as a -4 . Thus the fraction of energy density in GW increases. The earlier a mode re- 
enters the Hubble radius the more important this relative growth is, so the resulting inflationary GW spectrum today 
increases at high frequencies. This can significantly increase the amplitude of the inflationary GW in the frequency 
band of direct detection experiments, although one must check that the amplitude at higher frequencies satisfies 
the BBN bound (see section II Ap . Depending on the inflationary model, on the equation of state during the non- 
standard epoch, on the duration of this epoch and on when it occurs, this effect can be observ able by interferometric 
experiments, preferentially ground-based experiments since they operate at higher frequencies [l!4j |. It might also be 
observable by eLISA, e.g. if the GW at higher frequencies are sufficiently diluted by a matter-dominated stage taking 
place before the non-standard era with w > 1/3. Of course, if GW are produced by any other source between inflation 
and the end of the non-standard era, their amplitude would be amplified too and they would also be constrained by 
BBN. 

The predictions for the inflationary GW may also be modified if, in addition to the parametric amplification 
of quantum fluctuations, GW are also emitted classically by a non-zero anisotropic stress during inflation. GW 
produced inside the Hubble radius during inflation are first diluted by the exponential expansion of the universe, 
but their amplitude remains frozen afterwards once they leave the Hubble radius. Therefore, in order to minimize 
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the amount of dilution, the typical scale of GW production by a non-zero anisotropic stress during inflation must be 
sufficiently close to the Hubble radius. 

One example that has been recently discussed is the emission of GW by particles produced during inflation [l!5l - lll9l | . 
As the inflaton rolls down its potential, it provides a time-dependent effective mass to fields coupled to it. If such a field 
becomes effectively massless during inflation, particles of this field can be produced efficiently in a non-perturbative 
way. If particles are produced in sudden way at a given moment of time during inflation, they emit GW during a short 
amount of time before being rapidly diluted by the quasi-exponential expansion of the universe. This adds a bump to 
the standard inflationary GW spectrum, which can b e loc ated at any frequency toda y dep ending on the moment of 
time when the particles are produced during inflation jl!7| . In the models studied in [ll7| . this effect is however too 
weak to be observable by eLISA or ground-based interferometers. On the other hand, if the inflaton has derivative 
couplings to another field, this field may remain light throughout inflation and its particles can be efficiently produced 
in a continuous way. This occurs in particular when the inflaton <p couples to a gauge field through an interaction 
term of the form (fiF^F^, where F^ u is the field strength and its dual. Such an interaction is natural in models 
where the inflaton field is an axion. In that case, the continuous producti on of particles of the gauge field during 
inflation can significantly enhance the inflationary GW spect rum jf 17l . Illdj . However, these particles also lead to a 
highly non-gaussian contribution to the scalar perturbations [l!5j . so their production is constrained by the bounds 
on non-ga ussianities from observations of the CMB and the Large Scale Structures (LSS). In the models studied in 
[ll^ lTT9l]. the GW emitted by the gauge field may nevertheless be observable by ground-based interferometers, but 
not by eLISA, which probes scales that are closer to the ones probed by the CMB and LSS. 

Another example of inflationary GW produced by a non-zero anisotropic stress th at has been discussed in the 
literature occurs if a first-order phase transition takes place during inflation [l20l Il2ll ] . Again, one must check that 
the phase transition does not spoil the successful pre dicti ons of inflation. Different sources of GW from an inflationary 
first-order phase transition have b een c onsidered in [l2l| , where it was concluded that the signal is dominated by the 
collision of vacuum bubbles. Ref. |f 2f| also finds that the phase transition must complete relatively quickly, so that 
the parameter j3 introduced in section III Al is significantly larger than the Hubble rate H during inflation. The GW 
wavelength at the time of production, which is set by is therefore well inside the Hubble radius. As discussed 
above, the GW are then significantly diluted by the inflationary expansion before they leave the Hubble radius. A 
further suppression of the GW amplitude comes from the f act th at the energy density in the bubble walls is much 
smaller than the total energy density driving inflation. Ref. jl2l| then concludes that a negligible amount of GW is 
produced if the number of phase transitions is not very large. On the other hand, it is found in that paper that in 
inflationary models involving many first-order phase transitions, the GW could be marginally observable by LISA, 
which can perhaps be extended to the case of eLISA as well. Note that, because of the peculiar redshift of the 
modes during inflation, the shape of the GW spectrum differs from the one generated by first-order phase transitions 
occurring after inflation. 

GW produced by the parametric amplification of quantum fluctuations, but with an enhanced amplitude at high fre- 
quencies, have also been discussed in the context of the so-called pre-big-bang [l22j | and ekpyrotic/cyclic universe 123 | 
scenarios, which have been proposed as alternatives to the standard inflationary paradigm. Although qualitatively 
different, these two kinds of scenarios usually involve (in the Einstein frame) an epoch of accelerated contraction 
of the universe before the expanding era. Throughout the contracting phase, during which the tensor modes leave 
the Hubble radius, the background energy density increases with time, as opposed to the slow decrease taking place 
during slow-roll inflation. As a consequence, the GW spectrum today increases with the frequency. The amplitude 
is then negligible at CMB scales and is constrained by BBN at high frequencies. Depending on the evolution that is 
assumed in between the contracting and expanding phases, the GW may be observable by direct detection experi- 
ments, see [1, [f| for details. In between these two phases, the effective description breaks down and the dynamics is 
unknown. A better understanding of this regime is required in order to obtain firm predictions for the GW spectrum, 
as well as to assess the viability of these scenarios as alternatives to inflation. 



B. Gravitational Waves from the End of Inflation 



We now move on to GW produced after inflation. In fact, another privileged epoch for GW production is the end 
of inflation itself, or just after it, when the potential energy density driving inflation is converted into the thermal 
bath of the Hot Big Bang in the course of reheating. 

In many inflationary models, reheating starts with an explosive and non-perturbative decay of the inflaton con- 
densate into large, non-thermal fluctuations of itself and other bosonic fields coupled to it - a process called preheat- 



ing [I24j . The subsequent dynamics is dominated by a highly non-linear and turbulent-like evolution towards thermal 
equilibrium. The large field fluc tuations a mplified by preheating source a GW background that has been studied in 
several papers, see in particular |l25l4l31 |. Similarly to Eq. ([7]), preheating emits GW with a characteristic physical 
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wave- number fc* = ff*/e* at the time of production, where e* < 1 depends on the preheating mechanism and -ff* is 
the Hubble rate during preheating, which is directly related to the energy scale E- ln Q during inflation 13 . One then 
finds that the inflationary energy scale must be smaller than about E- ln ^ ~ 10 11 GeV for the peak of the spectrum to 
fall into the frequency range of ground-based interferometers (/ ~ 10 2 Hz), and smaller than E- m n ~ 10 7 GeV for the 
peak to fall into the eLISA frequency band (/ ~ 10 -2 Hz). Thus GW from preheating could be observable if inflation 
occurs at sufficiently small energy scales, which is complementary to the case of the GW generated during inflation 
itself. Note that an inflationary phase at low energies (possibly following a first stage of inflation at higher ene rgies) 
may have several cosmological advantages, e.g. preheating- induc ed baryo-g enesis |133j and magneto-genesis [134j, and 
the avoidance of dangerous relics such as moduli and gravitinos 

G W fr om preheating have been extensively studied with numerical lattice simulations as well as analytical methods, 
see |l28j for details. Similarly to Eq. ((5J), the peak amplitude depends essentially on e* and on the fraction of the 
total energy density that contributes to GW production. These two quantities can both be close to unity, which 
then maximizes the amount of GW produced. The peak frequency and the shape of the spectrum depend on the 
detail of the mechanism that is responsible for preheating, which in turn depends on the particular inflationary model 
that is considered. They also depend on the nature of t he fie lds that are produced during preheating, in particular 
gauge fields lead to specific features in the GW spectrum [130j . Two main classes of models have been studied so far: 
preheating after chaotic inflation and preheating after hybrid inflation. Another important cl ass o f models that may 
be interesting from the perspective of G W observation is preheating after small-field inflation jl37l | , but the resulting 
GW spectrum has not been studied with lattice simulations yet. In the following, we discuss in more detail the 
example of the GW from preheating after hybrid inflation. Chaotic inflation occurs typically at high energy scales, 
in which case the GW from preheating have a present-day frequency that is too high to be observable. 

In hybrid inflation models, the scalar field potential is of the form 

V=\ (\x\ 2 -v 2 ) 2 + 9 Y\x\ 2 <f> 2 + V s M) (66) 

where tf> is the inflaton and \ is another scalar field coupled to it. These fields can be complex or have an arbitrary 
number of components, but this has little effect on the GW produced by preheating itself. For > (j) c = y/Xv/g, 
where (j) c = \fX~vjg is called the critical point, the potential has a valley at % = 0. Inflation occurs when the 
inflaton slowly decreases in this valley due to the uplifting slow- roll term V sr (</)). When <j> = <fi c , the curvature of the 
potential becomes negative and the fields quickly roll towards the true minimum at <p = and \x\ = v. During this 
rolling, quantum fluctu ations of the fields have a negative effective mass squared and are exponentially amplified by 
a tachyonic instability [l38l 1 1 39j | . This quickly converts the homogeneous field energy into large fluctuations of the 
fields, which in t urn emit GW. Three different regimes of preheating and GW production in the model (|66[) were 
identified in |129| , depending on the couplings A and g 2 and on the inflaton velocity at the critical poin t (the latter 
depends on the particular inflationary potential V sr (4>))- For instance, in the regime g 2 <C A, Ref. [l29j | obtains the 
following estimates for the peak frequency and amplitude of the GW spectrum today 

f /peak ~ ^10 10 HZ 

For g 2 < A : { (67) 

\ h 2 Q qw \ . ~ io~ 5 Agv 2 

V y w \peak g z 

where in that case the condition Gv 2 < g 2 /X must be satisfied. 

GW from preheating in the model (|66p are observable in different regions of the parameter space, but this requires 
very small values of the coupling constants: typically g ~ vA < 10 -5 or A ~ 1 and g < 10~ 7 , se e Il29l for details. 
References to explicit pa rticl e-physics models of hybrid inflation wer e th is may occur are given in [129| . An example 
is provided by a variant [l4pf of brane inflation in a warped throat [l4ll ] , which is a prototype of hybrid inflation in 



string theory. According to [142j , the end of inflation in this class of models can be described by the effective potential 
(|66p with A ~ 1 and g ~ e~ Ab , where e~ Ab -C 1 is the warp factor at the bottom of the throat. In that case, a very 
small value of g may emerge naturally from the exponential warping of the throat, as small as g ~ e~ Ab ~ 10~ 16 if the 



Eq. is not directly applicable to preheating since by definition the universe is not in a thermal state and there is no definite 
temperature T* . In models where the full reheating process completes in less than one Hubble time, the reheat temperature is of the 
order of the energy scale at the end of inflation. Even if the full reheating process takes a much longer time to complete, the equation 
of state of the universe jumps usually quickly towards w as 1/3 at the beginning of preheating [1321 ] . In both cases, Eq. J7) with T* 
replaced by the energy scale at the end of inflation gives the correct order of magnitude for the peak frequency of the GW spectrum 
today. 
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hierarchy problem of the Standard Model is addressed a la Randall- Sundrum jl43| . The original model of |l4l| uses 
only a moderate warp ing, but it also suffers from a conformal-coupling problem. A model that avoids this problem 
has been proposed in [l4dl |. where e~ Ab ~ 10~ 16 — 10~ 13 and v ~ 1 — 1000 TeV. In such a case, Eqs. (|67| indicates 
that the GW spectrum from preheating can be observable by eLISA. 

In other inflationary models, inflation may also end with a vacuum (as opposed to thermal) first-order phase 
transition. In that case, bubble collision is an efficient source of GW, as already pointed out in Again, the 

peak frequency of these GW depends on the energy scale of inflation and can fall into an observable range if inflation 
occurs at sufficiently small energies. Recently, Ref. [l44| studied the GW produced by bubble collision at the e nd o f 
thermal inflation - an explicit inflationary model at low energy that may solve the cosmological moduli problem jl36| . 
In this specific model, inflation is followed by a rather long matter-dominated stage before reheating completes and 
the radiation-dominated era starts. GW from bubble collision are then diluted by this matter-dominated stage and 
their peak amplitude today is found to be /i 2 f2 gw ~ 10~ 18 in [l44j . which is too small to be observable by eLISA or 
ground-based interferometers. 



C. Gravitational Waves from the Hot Big Bang 



As discussed above, the present-day frequency of the GW produced by the non-perturbative decay of the inflaton 
condensate (i.e. preheating) depends on the energy scale of inflation. However, the inflaton is not the only scalar 
field condensate that can decay in such a violent and highly inhomogeneous way. Another example is provided by 
scalar fields that acquire a large amplitude along flat directions of the potential in super-symmetric theories. Th e 
GW produced by the non-perturbative decay of super-symmetric flat direction condensates were studied in |l45j . 
where it was found that they can be observable by ground-based interferometers. Indeed these condensates start to 
oscillate and decay when the Hubble rate becomes of the order of their soft super-symmetry breaking mass m. For 
m ~ TeV, this corresponds to a temperature T» ~ 10 10 GeV and Eq. ([7]) shows that the resul ting GW fall naturally 



in the frequency range of ground-based interferometers (e* is close to unity in these models 145j ). One difference 
with respect to preheating is that, while the inflaton condensate dominates the total energy density before it decays, 
flat direction condensates are typically sub-dominant. This may significantly reduce the GW amplitude, see the last 
factor in the RHS of Eq. ^/ . The GW are then observable for flat direction condensates that have a sufficiently large 
initial amplitude, which occurs for directions that are indeed sufficiently flat. 

Condensates along super-symmetric flat directi ons may a lso fragment into non-topological solitons, called Q-balls. 
The GW produced by this process were studied in [1461 Il47j | . A significant amount of GW can be produced for Q-balls 
with a large conserved charge. How ever, the decay rate of such Q-balls is small and they may come to dominate the 
energy density of the universe [l47j | . This leads to a matter- domin ated stage that further dilutes the GW produced 
from Q-ball formation. Taking this effect into account, Ref. jl47l | obtains GW amplitudes that are too weak to be 
observable by eLISA or ground-based interferometers. 



GW can also be produced by unstable domain walls 148Hl51 |. Domain walls originate from the spontaneous 



breaking of discrete symmetries. Once produced, they evolve towards a scaling regime, similarly to cosmic strings. 
However, stable domain walls would come to dominate the expansion of the universe and lead to cosmological disasters. 
This may be avoided if the discrete symmetry is only approximate. An example is provided by a real scalar field <p 
with low-energy potential 

v = \ (0 2 - v 2 ) 2 + e ^ (y-- 2 ) ■ (68) 

where e is a numerical coefficient called "bias" and the second term in the potential lifts the degeneracy between the 
two minima at <f> — ±v. This generates a pressure force on the walls that makes them collapse. The GW background 
from collapsing domain walls in the model (|68|) has been studied with lattice simulations in jl48l Il50j |. The GW 
spectrum can cover a very wide range of scales, from the width of the domain walls to the Hubble radius when they 
decay. This wide range of scales cannot be probed in lattice simulations, which must then be supplemented with 
analy t i cal es t imat es . 

Ref. jl49l | estimates the amount of GW produced by the collapse of domain walls in the following way. In a 
scaling regime, one expects the curvature radius of the walls and the average distance between them to be set by the 
Hubble radius, with the domain wall energy density evolving as pow ~ °" H where a ~ ^/Xv 3 is the domain wall 
tension. Stable domain walls would then dominate the total energy density at the time tdom when the Hubble rate 
is Hdom ~ G a. On the other hand, the bias term in Eq. (]68[) makes the domain wall collapse when the bias energy 
density V c ~ ei> 4 becomes of the order of Pdw-, which occurs at the time td ec when the Hubble rate is Hd ec ~ V e /a. 
The condition Hdec > Hd om is then required in order for the domain walls to collapse before they come to dominate 
the total energy density. On dimensional grounds, one expects that, in one Hubble time around tdec, the collapse of the 
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domain walls emits GW with physical wave- numbers £;* ~ H dec and energy density Ap gw ~ Ga 2 , which corresponds 
to the fraction Ap gw /p tot ~ G a 2 / (H^ EC / G) ~ (H dom /H dec ) 2 of the total energy density at the time of production. 
If the walls collapse during the radiation era, the amplitude of the GW spectrum today is then roughly estimated as 



( 2 

dom 



h 2 n gw (f dec ) ~ 1(T 5 [j^j (69) 

where the present-day frequency f dec is given by 

Uec~ (§f) 10- 4 Hz, (70) 

see Eq. 0. These estimates agree qualitatively with the numerical results of [l50j |. where it is also found that the 
GW spectrum is approximately flat at frequencies / > f dec . This suggests that GW from unstable domain walls could 
be observable by eLISA for T dec <; 10 3 TeV if H dom /H dec is not too small. The results depend strongly on the bias 
term in the potential, which in turn dep ends on the underlying particle-physics mo del. This has been studied in the 
context of gaugino condensation in [l49j and in the context of thermal i nflat i on i n jl5l| . 

Another source of GW that has been considered in several papers is the self-ordering of massless (or 



light) scalar degrees of freedom that may follow a global symmetry-breaking phase transition in the early universe. 
When a global symmetry is spontaneously broken, the massless Nambu-Goldstone bosons may lead to the production 
of GW at large, Hubble scales. An example is provided by a V-component scalar field <f> with O(N) symmetry and 
low-energy potential 



A 



V='-{\W-v 2 ). (71) 

During the phase transition, the scalar field settles in the true vacuum \<j)\ = v, but with different values of </> in 
different Hubble patches because the process cannot be correlated on distances larger than the horizon distance at 
that time. As the Hubble radius grows during the subsequent cosmological evolution, the scalar field relaxes into a 
uniform configuration over the expanding Hubble patches. This self-ordering can lead to a continuous production of 
GW at the Hubble scale from the anisotropic stress generated by the gradients of the Nambu-Goldstone bosons. This 
process shares similarities with the continuous "production" of GW when the tensor modes generated by inflation 
re-enter the Hubble radius. The GW produced by the self-ordering in the radiation-dominated era lead to a flat 
spectrum today with /i 2 f2 gw oc G 2 v A . Note that, compared to the case of inflation, the inflationary energy scale i?i n fl 
is replaced by the vacuum expectation value (VEV) v of the scalar field. 

Because the spectrum is proportional to G 2 v 4 , a very large VEV is r equi r ed fo r these GW to be observable. The 
GW produced by this mechanism ha ye be en calculated analytically in |l53l . Il54l | for the model ([7T]) in the large N 
limit, and with lattice simulations in jl55| |. Their results indicate that these GW could be observable by eLISA only 
if v > 10 17 GeV, which is larger than the maximal energy scale of inflation allowed by the current CMB bound. 
Nevertheless, the phase transition initiating the process could still occur after inflation if the potential is sufficiently 
flat (A -C 1 in Eq. (jTTj) ). On the other hand, as in the case of inflation, the spectrum can cover a very wide range of 
scales. If the source is active until GW at CMB scales are produced, then the spectrum can be subject to the COBE 
bound (see section O}. This can be avoided for instance if the relevant scalar modes are light instead of massless (e.g. 
pseudo-Nambu-Goldstone bosons if the broken global symmetry is only approximate). In such a case, we can expect 
that the source stops being active when the Hubble scale becomes of the order of the mass m of the light modes. This 
would introduce a low-frequency cutoff in the GW spectrum today, at a frequency corresponding to the wave-number 
fc~mat the time of production. 

Finally, GW are also produced at second o rder in cosmological perturbation theory by scalar modes re-entering 
the Hubble radius after inflation, see [156 158] and references therein. Since the scalar pertuibations at CMB scales 



are observed, the GW they produce are "guaranteed". These GW have however a small amplitude since they are 
produced at second order by small scalar perturbations. On the other hand, the scalar perturbations can have a 
larger amplitude at smaller scales, where they are much less constrained. This would increase the amplitude of the 
GW spectrum they produce, which in turn can provide complementary constraints on the scalar perturbations and 
on inflationary mod els. The GW prod uced by large density perturbations leading to primordial black holes (PBH) 
was then studied in [l59j | , see also |l60l | . It is found in particular that pulsar timing observations can probe PBH with 
~ 10 2 solar masses (which are candidates for intermediate-mass black holes), while eLISA could probe PBH in the 
mass range 10 22 — 10 25 g (which can play the role of dark matter). 
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V. CONCLUSION 

In this paper, we have investigated the scientific potential of eLISA in the area of cosmological backgrounds. 
We have considered several potential sources of gravitational waves that could fall in the frequency band of the 
instrument, focusing on the ones that are best understood: first-order phase transitions and cosmic strings. We 
have also reviewed several other sources that operate at different epochs in the cosmological history. In most cases, 
the frequency dependence of the gravitational wave spectrum can be estimated reliably, which allows to disentangle 
different cosmological signals from each other and from astrophysical or instrumental backgrounds. 

In Section |TT1 we analyzed the GW spectrum generated by first-order phase transitions, through the collision of 
bubbles and MHD turbulence. The analysis combines in a consistent way the most recent results of the literature, and 
provides therefore an up-to-date and realistic estimate of the GW signal. The analytic formulas of the GW spectra 
from bubble collisions and MHD turbulence, as a function of frequency and of the PT parameters a, /3/H*, r\ and T*, 
are taken from the latest analyses of Huber and Konstandin 2008 [26| and Caprini et al. 2009 [ll| . For the first time 
the bubble wall velocity Vb and the efficiency factor k have been evaluated beyond the Jouguet detonation hypothesis, 
adopting the complete model for bubble propagation developed in Espinosa et al. 2010 [23]. This causes the GW signal 
to depend also on the friction 77, acting on the bubble wall because of its interactions with the surrounding plasma. 
Moreover, the efficiency factor k has been constructed to account for both the gradient energy of the expanding wall 
and the bulk kinetic energy of the plasma in a continuous way, and it is dominated by one or the other depending on 
the strength of the PT: therefore, the contribution to the GW signal from bubble collisions and/or MHD turbulence 
is modeled consistently for each value of a. We also point out that, in a realistic model of the PT, the parameters a 
and cannot be considered as independent (see for example Huber and Kostandin 2008 We account for this 

fact when we evaluate the GW signal due to specific models of the PT, and show that it can decrease considerably 
the chances for detection. 

From our analysis of the GW spectra we provide detection forecasts in the case of a moderately strong (Fig. 2]) and 
of a very strong (Fig. [5]) first order PT. A key result is that eLISA will probe the presence of strongly first order phase 
transitions in the energy range from 100 GeV to 10 TeV or more, therefore complementing and even superseding the 
searches at colliders such as the LHC. We have given explicit examples, such as the holographic phase transition, 
where the signal lies well above the sensitivity of eLISA. There is of course no guaranteed detection but any such 
signal would point towards a strongly first-order phase transition: the frequency dependence is well understood, with 
a / 3 behavior at low frequency (due to the causal nature of process generating the GW), and a high-frequency tail 
given by the combination of the J™ 1 slope from bubble collisions and of the /~ 5 / 3 slope from MHD turbulence. If 
the phase transition occurs in the 10 to 100 TeV temperature range, eLISA might be the only way to identify it, and 
with a large signal-to- noise ratio if the transition is slow enough (lasting for a tenth of the Hubble time or more). 

In Section UHl we focused on the gravitational wave spectrum produced by a network of cosmic string loops. We used 
an improved model for the cosmological evolution and we studied the dependence of the results on the thermal history 
of the early universe. We checked that the gravitational wave background does not depend much on the spectrum of 
emission by individual loops (see Figs. [6l and [TO]) . We also found that removing the rare bursts is numerically relevant 
only for cuspy loops with a small initial size and in the region around the peak. As in the case of first-order phase 
transitions, the spectral dependance of the signal from cosmic strings is well determined. The high-frequency part of 
the spectrum (which is produced during the radiation era) is always almost flat over many decades of frequency (with 
small deviations from a perfectly flat spectrum arising from the variation of the number of relativistic species during 
the expansion of the universe). In the case of small initial loop sizes (a < TG/i), there is a well-distinct peak at low 
frequencies, with a very steep infrared tail and a tail at higher frequencies that depends (relatively weakly) on the 
spectrum of emission by individual loops. In the case of large initial loop sizes (a ~ 0.1), the peak is less pronounced 
and its infrared tail varies more slowly, as fl gw oc / 3 ^ 2 . 

An important feature of the background from cosmic strings is that it extends over many decades of frequency, 
because the source remains active during a long period of time. This may enable to detect the spectrum simultaneously 
with different experiments (pulsars timing measurements, eLISA, ground-based interferometers), which would improve 
the characterization of the signal and the ability to distinguish it from other stochastic backgrounds. We addressed 
the question of the detectability of the cosmic string signal by different experiments and in particular by eLISA. As 
shown in Figs. [9] and [T^l eLISA will be able to probe a wide portion of the parameter space. In a large part of the 
parameter space, the signal should be simultaneously accessible to at least one other experiment. In other regions, 
eLISA will be the only instrument capable of detecting the signal, thus opening a new window on the parameter space 
of cosmic strings. In the case of small loops, eLISA should be able to reach string tensions down to G[i ~ few x 10 -10 
for p = 1 and Gfi ~ few x 10~ 13 for p = 10~ 3 , in an appropriate range of loop sizes. For large loops, eLISA will be 
sensitive to G/j, ~ 10~ 13 for p = 1, and even G/i ~ 10~ 17 for p = 10~ 3 . 

Our main conclusion is that the scientific potential of eLISA regarding the detection of cosmological backgrounds 
is not significantly decreased from the original LISA mission. Indeed, the decrease in sensitivity is much lower than 
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the uncertainties that remain on the magnitude of the overall signal, due mainly to model dependence. This is why 
we have been careful in specifying the models chosen and the hypotheses made. As discussed in section ITV1 eLISA 
has also the potential to probe a rich variety of other physical phenomena and new ideas. 
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